Lowest-Landau-level theory of the quantum Hall effect: 
the Fermi-liquid-like state 



OO 
On 
OS 

u . 
< 



-a 
c 

o 
o 



> 
o 

OO 

-i— > 



I 

C 

o 
o 



X 



N. Read 

Departments of Physics and Applied Physics, Yale University, 
P.O. Box 208120, New Haven, Connecticut 06520 
(February 1, 2008) 

A theory for a Fermi-liquid-like state in a system of charged bosons at filling factor one is de- 
veloped, working in the lowest Landau level. The approach is based on a representation of the 
problem as fermions with a system of constraints, introduced by Pasquier and Haldane (unpub- 
lished). This makes the system a gauge theory with gauge algebra Woo. The low-energy theory 
is analyzed based on Hartree-Fock and a corresponding conserving approximation. This is shown 
to be equivalent to introducing a gauge field, which at long wavelengths gives an infinite-coupling 
U(l) gauge theory, without a Chern-Simons term. The system is compressible, and the Fermi-liquid 
properties are similar, but not identical, to those in the previous U(l) Chern-Simons fermion theory. 
The fermions in the theory are effectively neutral but carry a dipole moment. The density-density 
response, longitudinal conductivity, and the current density are considered explicitly. 



I. INTRODUCTION 

The so-called composite-particle view of the liquid 
states of electrons (or other charged particles) in two- 
dimensions in a high magnetic field |jj has been devel- 
oped gradually over more than a decade |^-p^]. Girvin 
Q proposed to develop a Ginzburg-Landau theory of the 
fractional quantum Hall effect, with an action for a com- 
plex scalar (boson) field and containing a Chern-Simons 
(CS) term to enforce the condition that the quantized 
vortices carry fractional charge. Girvin and MacDonald 
H introduced a singular gauge transformation and exhib- 
ited algebraic long-range order in a bosonic field. This 
transformation, which attaches <5-function flux tubes to 
particles (via a CS term in the action of the field theory) 
and so in general changes the statistics of the particles 
as in the theory of anyons ||l3|] , was then used in several 
theories, in conjunction with the mean-field approxima- 
tion of replacing the gauge field strength by its expec- 
tation value, to obtain a system in a different magnetic 
field. Thus, anyon superconductivity was discovered by 
mapping anyons in zero magnetic field to fermions filling 
Landau levels in a magnetic field H ; the Laughlin states 
jwj were described by mapping fermions to bosons in 
zero net magnetic field and then Bose-condensing them 
P|; the Laughlin and hierarchy 
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terpeted by mapping fermions to fermions in a reduced 
magnetic field and then filling Landau levels Hf7); the 
hierarchy states and the anyon superconductors in zero 
magnetic field were rcdcscribcd by hierarchical extension 
of the mapping to bosons, using duality methods M. At 
the same time, a lowest Landau level (LLL) treatment of 
the Ginzburg-Landau idea was developed ||[l0| , without 
using (^-function flux tubes, by attaching vortices to elec- 
trons to convert them to bosons; in this case, the bosons 
condense and have true long-range order. 



It has to be admitted that these ways of viewing the 
fractional quantum Hall effect produced little in the way 
of distinctive experimental predictions or explanations 
that were not already known by other methods, though 
interesting speculations concerning the phase transitions 
between the quantized Hall plateaus JlTj may be an ex- 
ception. The situation changed, however, following the 
discovery of an anomaly in the surface acoustic wave 
propagation at filling factor v — 1/2 (and less strongly 
at other filling factors, such as 1/4, 3/ 2) [ Q . This result 
speeded the development of a theory Q (to be referred 
to as HLR) for a case not included in the above list, 
in which fermions (electrons) are mapped to fermions at 
zero magnetic field and form a Fermi sea. In the simplest 
cases, this occurs for filling factor v = 1/2, 1/4, 1/6, 
The Fermi sea was predicted to be a compressible state 
that does not produce a Hall plateau, and the experimen- 
tal result of a longitudinal conductivity increasing lin- 
early with wavevector yl| was explained |ll[] . The Fermi 
surface, at which the fermions exist as genuine low-energy 
excitations, was observed through geometric resonance 
effects at v close to 1/2 in further surface acoustic wave 
experiments |?0| (as predicted explicitly in Ref . |ll| ) , and 
in other experiments [^1 22 1. (We should point out that 
for other filling factors in the fermion description, the 
fermions are dressed to become the fractionally-charged, 
fractional-statistics quasiparticles ]l4]- ^6| , p3[ , so are not 
observed as fermions.) 

In this paper, we return to the basic theory of the 
Fermi-liquid-like state. Recent work p4]-p6| has raised 
the possibility of changes in the way we think about 
the theory of the low energy excitations near the Fermi 
surface. In particular, these authors find contraints not 
mentioned in any earlier papers known to the present au- 
thor. At the same time, we may be motivated by trying 
to avoid the seemingly artificial CS approach, which be- 
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gins with a singular gauge transformation. Ultimately, 
it would aid our understanding to have more intuition 
about what drives the formation of the Fermi-liquid-like 
(and other) states. There are no flux tubes attached to 
the particles in reality; the background magnetic field 
remains essentially uniform in these states of matter. 
The approach begun in Ref. |l0| was intended to head 
in this direction. It uses LLL states only, so is valid 
in the (not entirely realistic) limit of interactions weak 
compared with u) c , and binds vortices to the electrons to 
lower the energy, thus forming the composite particles. 
Several implications of this approach were pointed out in 
Rcf. jl2) for the Fermi sea and the Bose condensate. 

The approach taken in the present paper avoids the 
CS approach. While it is perhaps not as simple-minded 
as one would want, it does make close contact with the 
work just cited fl^| . Here we start from an approach of 
Pasquier and Haldane (PH) |^,|^|, that gives an exact 
representation of the LLL problem in the case of charged 
bosons in a magnetic field at v = 1, where a Fermi liq- 
uid (FL) state is possible. Although our paper is long 
and fairly detailed, we can give a succinct summary of 
our results. The low-energy, long-wavelength theory is 
a FL coupled to a gauge field (not to be confused with 
the physical electromagnetic field). In contrast to the 
scenario arising ]Tl[ | in the CS (singular gauge transfor- 
mation) approach, there is no CS term in this low en- 
ergy theory. Consequently, the gauge field is said to be 
"strongly-coupled" and one of its effects is to enforce con- 
straints that agree with those of ]24|-|2"6|] . This in turn 
has the effect of making the fermions uncharged, but 
they are left with a subleading coupling to electromag- 
netic fields through a dipole moment. The interplay of 
this moment with the transverse part of the gauge field 
leads to a finite compressibility, in spite of the neutral- 
ity of the particles. It also leads to the CS equations, 
that relate the curl of the vector potential to the den- 
sity, and a similar equation for the current, still being 
valid, in spite of the absence of a CS term in the action, 
in agreement with Ref. [O. In general, the good agree- 
ment with experimentally-observed phenomena achieved 
in the theory of HLR is not spoilt in the present theory. 
Nonetheless, the detailed structure of this FL-like theory 
is modified. While the theory is developed here for v = 1 
bosons, there are many indications that the results are 
more general. These include the derivation in Ref. p4| ] 
for general number of attached flux. 

Sec. [n] contains a more detailed review of previous 
work, and a more detailed overview of the paper. In Sec. 
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we explain the formalism due to Pasquier and Hal- 
dane that will be used in this paper. In Sec. [Tv| , we per- 
form explicit calculations of response functions, including 
those for the constraint operators, and interpret the re- 
sults in terms of a strongly-coupled gauge field. In Sec. [V| 
we outline the extension of the results to all orders, and 



sion. Appendix |A| discusses some details of the formal- 
ism, including the noncommutative Fourier transform, 
and Appendix [B] indicates how a Hubbard-Stratonovich 
transformation can be used. 



II. REVIEW AND OVERVIEW 

In this Section, we review some of the background nec- 
essary for the discussion in this paper. We begin with the 
U(l) Chern-Simons (CS) fermion approach developed in 
Rcf. jll). The Fermi liquid-like state proposed in that 
paper is the main topic of the present work; however, 
we will not review the relation to experiments. In Sub- 
sec. 



[I B we review "physical" pictures which are based 



on consideration of the wavefunctions of the system, as 
opposed to field theoretic methods. 
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In Subsec. 

review recent work which attempts to push the U(l) CS 
approach down to a low-energy effective theory in the 
lowest Landau level (LLL). Finally, in Subsec. II D, we 



give a brief overview of the main results and of the layout 
of the remainder of the paper. 



A. U(l) Chern-Simons fermion theory 

In this approach the particles arc represented as 
fermions with a (^-function of flux attached, of strength 
an integral number (f> of flux quanta <i>o. Then the under- 
lying particles must be bosons when cj) is an odd integer, 
and fermions when cf> is even (for noninteger <f>, the un- 
derlying particles must be anyons). We will reserve the 
term "particles" for these original particles, and refer to 
the transformed particles as "fermions" or "quasiparti- 
cles". The imaginary time action (see, e.g. \LM, to be 
referred to as HLR) is (in the gauge where V • a = 0) 



S = 



dr d 2 j 



1 



ip* (J?- - ia - f?j tp 



— |(-iV a AUI 2 ^a V A a 

2m 27T0 

+ dTd 2 rd 2 r'V(r-r')^{r)^(r')^(r')^(r). 

(2.1) 

Here i\> is the field operator for the fermions, rather than 
that for the underlying particles, which could be fermions 
(electrons) or bosons. We will use the notation (note the 
use of the summation convention for repeated Greek in- 
dices) 



a A b = e„ v aJ) v 



(2.2) 



provide some general discussion. Sec. VI is the conclu 



for a cross product of vectors a, b in two dimensions, 
/i, v, . . . = x, y to label the two components, and 
s xy — 1 for the two-dimensional alternating 
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tensor. We have set % = 1 and, starting with gaussian 
units, we have absorbed — e into the scalar potential and 
electric field, and (— e/c) into the vector potential and 
magnetic field, so the charge of the particles is one and 
the flux quantum is 2tt. The uniform background mag- 
netic field isVAA = £>>0 which corresponds to the 
negative z direction (in the three-dimensional sense) in 
conventional units. We choose the unit of length so that 
the magnetic length £g 2 = B = 1. It will also be conve- 
nient to write Aa for the vector whose components are 
(Aa) M = £ MI/ a„; then a Ab = a A b. 
Varying ao in the action leads to 



VAa = -27TC 



(2.3) 



where p(r) = ^(r)ip(r) is the number density both of the 
Chern-Simons fermions and of the underlying particles. 
When the filling factor v = 2np/ B is l/<fi (where p is the 
average density), there is no net field for the fermions, 
and within a mean-field approximation, a Fermi sea 
ground state is possible. 

The leading approximation for the linear-response 
functions is the random phase approximation (RPA), in 
both the gauge field ao, a and the Coulomb (or other) in- 
teraction V(r). In Fourier space the full density-density 
response function is then |fl| , before any approximation, 



Xpp 



1 + V(q) X 



ITT 
PP 



and in the RPA Xpp — X 1 " i where 



Xo = 



Xo 



1 - (27r£) 2 XoXo79 2 



(2.4) 



(2.5) 



Here x 1IT is the response function which is irreducible 
with respect to the interaction V only (i.e., diagrammat- 
ically, it does not become disconnected when a single 
interaction line is cut), while xo is the density-density 
response for the non-interacting sea of fermions of mass 
m (the bare or band mass), and Xo" is the transverse 
current-current response, of the same Fermi sea, includ- 
ing the constant "diamagnetic current" term. In the limit 
where first the frequency u> and then the wavevector q 
tend to zero, we have 



Xo = to/2tt, 

Xo o 2 /127rm, 



and hence 



On _ h 
— = km x p 



;m) = 



m/27r 
l + (^ 2 /6' 



(2.6) 
(2.7) 



(2.8) 



(For a long-range potential, i.e. one that is divergent as 
q — > 0, this is the appropriate definition of the compress- 
ibility dn/dp. For a short range interaction, one would 



use Xpp m place of X P ") Thus the theory predicts that 
the system is compressible. Note however that the ap- 
proach describes the properties that the system has if it 
is in the phase described. For a highly-correlated system 
such as particles in the lowest Landau level, it is difficult 
to find any approach that can accurately predict, for a 
given Hamiltonian, in which phase the system will be. 
For example, an alternative phase that is possible at the 
same filling factors as the Fermi liquid (FL) is the Pfaf- 
fian state |^8| , which is believed to be incompressible |29| . 
Nevertheless, the question of the properties of the Fermi 
liquid state — which has a Fermi surface in the excitation 
spectrum for the fermions — is well-defined. 

For the conductivity, the general statement |?0| is that 
the resistivity tensors add: 



P — Pes + Pi,, 



(2.9) 



where pcSp.v = 2TT(f>e^, coincides with the Hall resistiv- 
ity at v = l/<p and p^^y is the resistivity tensor of the 
fermions, the inverse of the conductivity tensor which is 
related to the current-current response function that is 
irreducible with respect to both the interaction and the 
gauge field. In the RPA, using the Drude approximation 
to include impurities, one has Jll| at q — > 0, then u) — > 0, 
Pi>p,u — S^v/a^xx where a^ xx is the usual Drude result 
for the Fermi sea in zero magnetic field with impurity 
scattering. There is also an unusual scattering mecha- 
nism |ll]|n| in which the fermions scatter off the static 
vector potential <5a induced in the Chern-Simons gauge 
field by a variation in the density of particles produced 
by the impurity potential, since V A <5a — —2'K(j)8p. 

The effects of interactions and gauge field fluctuations 
beyond RPA would be expected to have a variety of ef- 
fects. By analogy with the Landau-Silin treatment of 
fermions with a long range interaction, one would expect 
that when both the long-range part of the interaction (if 
any) and of the Chern-Simons gauge field are extracted, 
by considering responses irreducible with respect to both 
the interaction and the gauge field as above, the remain- 
ing effects can be handled to all orders by renormaliz- 
ing parameters, and the leading long-wavelength effects 
expressed in terms of Landau interaction parameters Fg, 
and an effective mass m* . Since the system is translation- 
ally and Galilean invariant (in the absence of impurities), 
the latter mass must satisfy the usual relation |il],p2]| 

m*/m = 1 + (2.10) 

(details of our two-dimensional normalization of the Lan- 
dau parameters such as F\ are given later). In addi- 
tion, in the limit where the cyclotron energy uj c — l/m is 
large compared with the typical interaction strength be- 
tween particles, V(p~ 1 / 2 ), (e.g. as to — > 0), the dynamics 
should be governed entirely by the interactions, and so 
1/to* should scale with the interaction strength, and be 
of order the typical interaction strength up to numerical 
factors. 
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This expectation that the theory would be a renormal- 
ized Fermi liquid, coupled to the long-range interaction 
and the gauge field, turned out to be too naive, however. 
The fluctuations of the gauge field have singular effects 
that appear to cause a partial breakdown of the Fermi 
liquid picture pT| |. The effects of such fluctuations were 
evaluated in leading order in the RPA gauge field prop- 
agator in HLR (the small parameter being cj>, with the 
background magnetic field being adjusted such that the 
net field seen by the fermions on average was zero, for 
any value of <fi, i.e. the filling factor was always 1/0; re- 
call that for generic values of <fi the particles are anyons). 
The main effects were, first, that the propagator itself 
shows the appearance of a mode at the cyclotron fre- 
quency l/m, which carries all of the f-sum rule spectral 
weight to order q 2 . Thus this mode is the physical cy- 
clotron mode. The virtual excitation of this mode, which 
is the longitudinal part of the gauge field, led, in first or- 
der, to a contribution to the fermion self energy that 
was logarithmically infrared divergent. The effect could 
plausibly be exponentiated to give for the quasiparticle 
residue Zp of a fermion at the Fermi wavevector kp, 

Z f ~ L-* /2 , (2.11) 

where L is the system size (or, presumably, \k — kp\^l 2 as 
k approaches kp for infinite L). This would correspond 
to the Girvin-MacDonald (GM) power law H, gener- 
alized to the fermion case; in particular, the exponent 
should be exact. This is supported by further analysis 
of these fluctuations which, similarly to the boson case 
p3|] , lead to a factor Y\.i<j \ z i ~ z j\ ■> times a gaussian, in 
the ground state wavefunction of the fermions (the result 
for the fermion case is widely known but does not appear 
to have been explicitly published). This in turn leads to 
the GM power r~^l 2 as a factor in the equal-time Green's 
function of the fermion [§4j 

(^(r)^(o)) ~ r -^l 2 +il 2 ) sin(fc F r - tt/4) (2.12) 

and correspondingly to the above result for Zp (see also 
Ref. HH). (The GM power law in the composite boson 
case has also been recovered field theoretically in [p6fl-) 
Related effects were also found in the work of Shankar 
and Murthy [M, to which we shall turn shortly. In the 
work of HLR and others, it was assumed that the vanish- 
ing quasiparticle residue for the original CS fermions was 
of little significance, since as with many similar effects in 
field theory, in particular the nonsingular quasiparticle 
residue in an ordinary Fermi liquid, it cancels in physi- 
cal response functions that measure quasiparticle prop- 
erties. However, the recent results to be reviewed below, 
and those of the present paper, suggest that things are 
not quite so simple, and rather than just ignoring these 
effects on the assumption that they cancel, the longitu- 
dinal mode should be integrated out "exactly" to obtain 



an effective field theory, before proceeding to the effects 
of the other lower-energy fluctuations, such as the trans- 
verse fluctuations. 

The fluctuations in the transverse part of the gauge 
field have received more attention (due to the CS term, 
there are also cross-terms that mix the longitudinal and 
transverse fluctuations; however these are assumed to 
have some intermediate significance). The first-order self 
energy contains power-law infrared-divergent terms for 
the case of a short range interaction, which are weakened 
by the presence of a long-range interaction because the 
latter suppresses density fluctuations which correspond 
to fluctuations of the transverse CS vector potential a. 
For the 1/r Coulomb interaction, the effects become log- 
arithmic, and for an interaction which is longer range 
than 1/r they become finite. In the Coulomb case, the 
structure of the effects is similar to those in an electron 
gas coupled to the transverse part of the ordinary elec- 
tromagnetic field (since there is no CS term in this case, 
these effects are not weakened by the Coulomb interac- 
tion, but are always logarithmic — however, they are ex- 
tremely weak in practise) |37],[38|]. In both of these sys- 
tems, it can be argued by treating the self energy self- 
consistently [ ^7|JTl|j39| ] that the effects lead to an effective 
mass diverging as m* ~ — In \k—kp\, a quasiparticle scat- 
tering rate ~ — \e% — In \e* k — (where e* k is the disper- 
sion relation that corresponds to the stated behavior of 
the effective mass near fc_p), and a quasiparticle residue 
Zp ~ — In |fc — kp\ (the latter would be in addition to the 
effect of the longitudinal fluctuations described above). 
These results suggest that while the effective mass di- 
verges at fci?, the quasiparticles remain just marginally 
well-defined due to the reciprocal logarithm in the decay 
rate, and thus the system is a "marginal Fermi liquid" . 
For longer-range interaction, there is no such breakdown 
of Landau Fermi liquid theory, and for the extreme case 
of V{r) ~ lnr, the scattering rate recovers its usual form 
~ (efc — /i) 2 (all these results are for zero temperature). 

There are many other studies of this [^p|-^6| , often with 
conflicting results. We believe that the correct results are 
those that agree with the above scenario of HLR for the 
behavior of the effective mass, etc. 

If we are not too concerned about the latter effects 
of transverse gauge field fluctuations, for example if we 
consider an interaction longer-range than Coulomb, or 
in the Coulomb case neglecting the logarithmic effects 
in view of how slowly they diverge at kp, then we are 
led to a physical picture of what to expect from the sys- 
tem to all orders in the fluctuations. It is essentially the 
Landau theory with due regard to the long-range effects, 
as described above, and thus retains the CS structure 
present in the RPA. For the density-density response, the 
responses xo and Xo that appeared in the RPA will there- 
fore be replaced by renormalized versions, and according 
to this scenario, we then expect that, in the limit that 
gives (for example) the compressibility, xo an d Xo" that 
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appeared in the RPA will be replaced by renormalized 
versions, m*/[2w(l + F )] and q 2 Xd respectively, where 
is a renormalized long-wavelength Landau diamag- 
netic susceptibility, which is a non-Fermi-liquid quantity 
as it involves derivatives at the Fermi surface. Explicitly, 
we expect 



dn 



m 



2tt(1 + F ) - (2^) 2 X *m* 



(2.13) 



(We expect that F diverges the same way as m* so that 
the renormalized version of Xo remains finite @|[§.) 
Thus the system remains compressible in this scenario. 



B. Physical pictures 

In this subsection, we review aspects discussed in Ref. 
[ 12 1 , which was in part an elucidation of Ref. |l(| (see also 
[47|). The approach is based on the wavefunctions of the 
particles, which are assumed from the beginning to be 
in the lowest Landau level (LLL) . To lower the repulsive 
interaction energy, each particle would like to bind to (f> 
vortices, which at v = l/4> leaves no vortices left. (Note 
that in the LLL, the number of zeroes in the wavefunc- 
tion of each particle is equal to the number of flux quanta 
threading the relevant area, and that a vortex means a 
simultaneous zero in the wavefunction of every particle 
other than the one under consideration.) For the same 
choices of statistics of the particles and of filling factor as 
before, the bound states behave as fermions in zero net 
magnetic field (this statement again involves the mean 
field assumption that the average density of particles is 
uniform, as we will see). 

To make the idea concrete, we may consider trial wave- 
functions in which the fermionic bound states occupy a 
Slater determinant of plane waves, or spherical harmonics 
on the sphere (these resemble Jain's trial wavefunc- 
tions |(|, except that the fermions are in zero effective 
magnetic field) 

*(^,... ) z JV )=P LLL detM y Ylizi-Zj)*. (2.14) 

i<j 



Here we write the wavefunction on the sphere |jl5|| , with 
Zi = 2Rvi/ui the complex coordinate of particle i in 
stereographic projection to the plane. Only the poly- 
nomial part of the wavefunction is shown, as indicated 
by the tilde on ^ . The full wavefunction is recovered by 
multiplication by JJ^l + |z l | 2 /2i? 2 ) _(Ar * +2)/2 , and this 
must be done before integration of the Zi coordinates 
over the complex plane, to give the correct integration 
measure, in particular when applying the LLL projec- 
tion operator "Plll • In the limit where the radius R and 
the number of flux quanta through the surface of the 
sphere go to infinity with the field strength fixed, the 



_i y u.i2 

non-polynomial factor approaches the usual e 4 ^< 1 1 . 
Mij are the spherical harmonics of angular momentum 
Li, Mi for the jth particle, or can be replaced by plane 
waves e iki ' rj in the plane. The Li, Mi (or k,;) can be 
chosen to fill the Fermi sea to obtain a trial ground state 
p8| . Different choices of the sets of Li, Mi do not give 
orthogonal states in general, except when the total angu- 
lar momenta differ. Note that apart from the projection 
to the LLL, the wavefunctions have the form that would 
be expected from the CS approach, on including the fluc- 
tuations at the RPA level that produce the amplitude of 
the Laughlin-Jastrow (LJ) factor in the wavefunction, as 
noted above. 

The fermionic bound states or "quasiparticles" de- 
scribed here are created by operators of the form iplU^, 
where ipl creates a particle in the LLL, and U(z) — 
Yli( z i — z) is Laughlin's quasihole operator |Q, which 
creates a vortex jL0|| . As for the wavefunctions, this dif- 
fers from the CS fermion operator ip^ by including the 
amplitude of the quasihole operator, and not just the 
phase (like the wavefunctions, it should also include a 
non-polynomial factor in z, which we have suppressed 
here). Consequently, like the corresponding boson oper- 
ator [jlO| , its equal-time Green's function is not expected 
to include the GM power-law factor r~^/ 2 ; this has been 
confirmed by calculation |Q. Since at v = l/</> the 4> 
vortices induce a hole in the density of the other par- 
ticles that contains a deficiency in the particle number 
of exactly unity, there has always been a temptation to 
say that the bound states formed this way are neutral 
objects. This should be contrasted with the CS fermions 
and bosons, which carry particle number unity. 

The plane-wave factors, in the flat space limit, can be 
rewritten using {l£|^0| (see also Appendix |A|) 



pLLLe lkr 'P L LL = e 



k R; 



■m 2 



(2.15) 



where R; is the guiding-center coordinate of particle i, 
which has no matrix elements between states in differ- 
ent Landau levels. The operator Kj = — A R^ is the 
pseudomomentum that generates magnetic translations 
of particle i. Thus, the plane- wave factors in the Slater 
determinant can be replaced by e lk Ri and each such fac- 
tor displaces the ith particle by Ak (in units where the 
magnetic length is one) from its vortices. This picture 
of particles bound to vortices but displaced by Ak from 
their center has several consequences J^j . 

The first consequence is that, if we consider the in- 
teraction of the particle with the vortices (or correlation 
hole) to which it is bound (neglecting the exchange effects 
due to the latter being constructed from other particles, 
indistinguishable from the first), then for k = 0, the par- 
ticle is precisely on the vortices as in the Laughlin states, 
and for k ^ it is displaced by Ak. Consequently, the 
energy should increase, and the interaction between the 
particle and its vortices becomes an effective kinetic (i.e. 
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k-dependent) energy for the fermion, which is the ori- 
gin of the effective mass at the Fermi wavevector, and 
scales inversely with V. A formula for this energy can 
be found for the analogous boson case in Ref. ]lO[j . No- 
tice that the displacements in the Fermi sea ground state 

are bounded above by kp = J 2/</>, which is much less 
than the typical distance between neighboring particles 

which is of order ~ Thus for <\> > order 1, which 

is the case of interest when the particles are bosons or 
fcrmions, not anyons, the displacements do not unduly 
perturb the bound states. 

Second, if we accept that the fermions are neutral, 
then their leading coupling to the electric potential is 
through a dipole moment Ak. It is important to realise 
that the wavevectors of the fermions contribute to the to- 
tal momentum of the system, which is a conserved quan- 
tity. One might imagine that the dipole moment could 
be renormalized by effects not yet included, or that the 
vortices might not all be at the same point as we have 
implicitly assumed. Indeed, when the underlying parti- 
cles are fermions, the wavefunction must have one vortex 
exactly on every particle, because of antisymmetry. This 
will not affect the dipole moment, because the plane- wave 
factors must produce the displacement shown, and when 
the particles are fermions, this is accomplished by dis- 
placing the other vortices further to compensate for the 
one that is not displaced at all. Also, if the vortices are 
viewed as point objects, then their relative displacements 
can only produce multipolc moments of even order, and 
not a contribution to the dipole moment, which is deter- 
mined by the displacement of the particle from the center 
of mass of the vortices. Thus the dipole moment is not 
renormalized. A more rigorous version of this argument 
will be given later in this paper. 

Third, when the <f> vortices are dragged around adia- 
batically, they pick up a Berry phase factor [p3| which can 
be interpreted as a vector and scalar potential governed 
by the particle number and number drift-current densi- 
ties, p and j |l^,[l2],[|7|] . This means that the fermionic 
bound states experience, in addition to the electromag- 
netic A and Aq, also a, ao given by 



V A a = -27r</>p, 
-a — Vao = 2-K<j) A j. 



(2.16) 
(2.17) 



These have the form of the equations in the CS fermion 
approach, but it is important to emphasise that they have 
been obtained |T(||l2|,[l7| without the use of 5-function 
fluxes attached to the particles, and that they still in- 
volve the physical density and current, which cannot be 
identified with the density and current of the fermions 
because the latter are (or may be) neutral. 

In Ref. Q, these were used as an alternative ap- 
proach that was stated to be equivalent to the CS ap- 
proach, and the neutrality of the quasiparticles was not 



invoked. It was felt that, although the fermions and 
bosons appear neutral, the situation might be like that 
in the usual electron gas problem with a Coulomb in- 
teraction, where at low energies the quasiparticles are 
neutral in their couplings to external longitudinal elec- 
tric fields because of screening, but in the Fermi liquid 
viewpoint, one nonetheless views the fermions as having 
charge unity, and the low-energy behavior of the Fermi 
liquid itself produces the screening effects, in the limit 
uj/q — > in the response functions. In the opposite limit 
uj/q — > oo, the charge of the quasiparticles does show up, 
in the conductivity (and also in the transverse response 
in both regimes). However, recent work to be discussed 
in the next subsection, and the work in the present paper, 
suggests that in the quantum Hall effect context, we can 
in fact obtain a consistent picture in which the quasi- 
particles have only dipolar couplings to external fields. 
The obvious question is then whether the Fermi liquid 
is still compressible. We will answer this question in the 
affirmative. 



C. Recent approaches to the LLL 

Several recent works have taken up the outstanding is- 
sues discussed in the previous subsections. They are con- 
cerned with obtaining results for the Fermi-liquid state 
including the effects of all the particles being in the low- 
est Landau level, or as would seem to be at least roughly 
equivalent, including the effects of the amplitude of the 
correlation factors produced by the zero-point fluctua- 
tions of the cyclotron-frequency longitudinal modes of 
the CS gauge field. The aim of such work is, of course, 
to test the validity of the results of HLR. Different ap- 
proaches have been used. Shankar and Murthy (SM) |24|1 
base their work on the U(l) CS fermion field theory ap- 
proach, however they work in a Hamiltonian formalism, 
and aim to eliminate the cyclotron variables by canonical 
transformation, rather than by resummation of perturba- 
tion theory. The cyclotron modes are represented as os- 
cillators whose zero-point motion produces the amplitude 
of the LJ factor in the ground-state wavefunction. How- 
ever, when fermions are excited to different k states, the 
oscillators must adjust to a displaced ground state, and 
this seems to reproduce many of the effects of the corre- 
lation hole discussed in the preceding Subsection, as well 
as other effects connected with the cyclotron mode and 
the projection to the lowest Landau level. D.-H. Lee |25| ] 
uses duality methods, which are good for representing 
vortices. In his approach, the particles are fermions at 
v = 1/2, but, in view of the single vortex exactly on each 
particle because of Fermi statistics (for LLL wavefunc- 
tions), they can be represented as bosons at v = 1. In 
these two works, only the leading long-wavelength effects 
can be treated. Pasquier and Haldane (PH) j2(| use a 
method that is valid only for </> = 1 (that is, the particles 



G 



are bosons at v — 1), and represents the LLL problem 
exactly, through equations valid for all wavelengths. A 
version of their method will be described in the next sec- 
tion and used extensively in this paper. 

All these groups arrive at the following points in com- 
mon. The LLL physics is described by Fermi fields c, 
with canonical anticommutation relations, and the phys- 
ical states must obey the operator constraints for each 
wavevector q, 



d 2 k 

W) 



2 L k- 



5 c k+i q (l- ^*Aq) 



+ 0{q 2 ) - p(27r) 2 <5(q) = 0. 



(2.18) 



In SM and Lee, the 0(q 2 ) terms are unknown, and in 
SM the constraints are further restricted to apply only 
for q less than a cutoff Q whch is chosen to equal hp. In 
PH, the terms higher-order in q are known. The physical 
particle number-density operator reduces to the form 



p(q) = p(2tt) 2 %) 



cPk 



*kAqc£_, c k+ i q , 



(2.19) 



again to leading order in q, on using the constraints. Note 
that this is the Fourier transform of a dipolar or polariza- 
tion expression for the density, p = p — V • P, where the 
polarization P is that due to a dipole moment of Ak on 
a fermion of wavevector k (this semiclassical way of de- 
scribing it will be quite useful; compare the discussion of 
fermions with a fairly well-defined wavevector and posi- 
tion in Fermi-liquid theory, which can be better described 
formally by the Wigner distribution function). Lee dif- 
fers from the other authors and from Ref. p2[ in finding 



an extra factor of 1/2 in the right-hand side of Eq. (2.19); 
the origin of this 1/2 is not clear to us. 

A result for the effective mass was obtained as follows. 
Beginning from the interaction Hamiltonian that is all 
that is left when the kinetic energy of the particles has 
been quenched, 



~<l 



(27T) 



V(q) : p(qM-q) : 



(2.20) 



where colons : . . . : represent normal ordering, the normal 
ordering is then dropped as it produces only a constant 
proportional to the number of part icles . The density is 
then replaced by the form in Eq. ( 2.19| ). When this is 
written in first quantization it becomes 



-ffint — 



2^ 



d 2 q 



V(g)q A kiq A kj 



(2.21) 



On taking the i = j term of this expression, they obtain 
an effective kinetic energy due to interactions, 



]>> 2 /(2to*) 



(2.22) 



where the effective mass is given by 
1 f d 2 q 



1/m* 



(27T} 



:V(q)q 2 



(2.23) 



which has the form of the dipole-moment-squared term 
in the self-interaction energy of a dipole, and if the q in- 
tegral is cutoff as in SM, the density profile is smeared as 
it would be in the correlation hole. It is therefore similar 
to the proposal of Ref. 0Jl|]. 

For the density-density response function, these au- 
thors find, using the dipolar form of the density, 



Xpp(q,0) = (p(q)p(-q)) 



(PP) = q 2 pm* + 0(« 4 ). 

(2.24) 



In the last step, the transverse momentum- momentum 
response function of the Fermi gas with effective mass 
to* was used. In these calculations, the constraints ( 2.18) ) 
were either ignored [p6[ , or were handled by introducing 
functional-integral representations of (5-functions of the 
constraints, which were then treated in the RPA [pi[ ; 
the results take the form stated in either case. 

If this last result is taken seriously, it implies that the 
system is incompressible. However, SM state some reser- 
vations about the calculation, because of the way the 
constraint was handled. They suggest that the symme- 
try of the Hamiltonian under translations of the wavevec- 
tors of all the particles could lead to cancellations and to 
factors of 1 jq 2 that could restore a finite compressibility 
to the system. This proposal is very close to the results 
that will be obtained in the present paper by a system- 
atic treatment of the constraints. While this paper was 
being completed, a short comment | ^T[ and a revised ver- 
sion of Ref. p5| ] appeared which use the same symmetry 
just mentioned and obtain results very similar to some 
of ours below, including the fact that the system is com- 
pressible. We will comment further on the relation of the 
symmetries being used in Sec. [v| 



D. Overview of the results of the present paper 

Here we describe results of the present paper. First 
we give a simple discussion of our central result, for the 
density-density response function. With the benefit of 
hindsight, using arguments that are justifed by the more 
detailed and formal calculations below, th e res ults can in 
fact be obtained from the results of Sec. 
describe the results of this paper. 



II B. Then we 



The dipolar form of the density in Sec. |II B| can be 
expressed as 



p(r) = p - V A ; 



(2.25) 
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where g(r) is the momentum density of the fermions, 
since P = Ag. On the other hand, we also have 



p(r) =p-pVA(a + A). 
This suggests that we write 

a + A = g/p 



(2.26) 



(2.27) 



in general, even though the above argument only implies 
this for the transverse part of a. This equation suggests 
there is a gauge- invariant current j , which is not the 
physical current, such that (for excitations near the Fermi 
surface) , 

f = |-^[c f Vc- (Vc'jc] - (a + A)c t c| /m* 

= [g(r)-(a+A)p fl ]/m*, (2.28) 

which is required to vanish, j R = 0. Assuming that the 
"density" p R — c^c is just p, this is equivalent to Eq. 



(2.27). Indeed, vanishing current would be consistent 
with such a constraint, p R — p, if they together obey a 
continuity equation, 



dp R /dt + V ■ f = 0. 



(2.29) 



This involves the longitudinal part of the current, so we 
have an argument for both parts of Eq. ( 2.27 ). The con- 
dition p R — p should of course be viewed as the long- 
wavelength version of the constraint found by SM, Lee, 
and PH. 

The gauge-invariant form of the "current" j R encour- 
ages us to consider an effective Hamiltonian 



H ^ = 77^-1 /d 2 r|HV-a-A) C | ; 
2m* J 



(2.30) 



which, apart from higher covariant derivatives of c, c\ 
contains no other terms in a, not even a CS term. Thus a 
is a strongly-coupled gauge field and varying H c g with re- 
spect to a yields j R = 0. Then, neglecting other terms in 
H c fj , we can use the RPA, or the following self-consistent 
field argument, to find the density-density response func- 
tion. From the form of the density, an external scalar po- 
tential couples to V Ag. The irreducible density response 
contains two parts. The first part is from the transverse 
momentum- momentum response function of the gas with 
mass to*; it is the part found by SM, Lee, and PH 24 2(J 
(Lee has since revised this result J25| ) . The second is the 
response of the same gas to the induced vector poten- 
tial a. (In both responses, the constant "diamagnetic 
current" term is absent.) Thus 



Xpp = (p + m *Xo )(<7 m * + iq$a±_) 



(2.31) 



where in the last factor the two terms arise from the 
two parts just described, and Sa± is the response in the 



transverse vector potential to the perturbation, and is 
therefore given by 



iqSa± = x p T p/p- 
From these self-consistent equations we find 



X, 



-p(p + m*Xo )q 2 /Xo, 



(2.32) 



(2.33) 



which is exactly the result we obtain in this paper. 
This yields for the compressibility dn/dp — —p 2 /Xd > 
0, where Xd 1S the diamagnetic susceptibility for this 
fermion gas. This result differs from that in the scenario 
based o n the U(l) CS approach, described at the end of 
Subsec. II A . Several other observables are similarly in 
close, but not always exact, agreement with the scenario 
based on HLR, described above. 

In this argument, we neglected the Landau parameters. 
These can be included without significantly changing the 
results. However, the Landau parameter F\ should not 
be added, since it is already included in the gauge field ef- 
fects. The strongly- coupled gauge field in the Fermi liquid 
is equivalent to a Landau parameter Fi — —1, provided 
to* > 0. Thus we are led to a scenario in which the Fermi- 
liquid-like state has many FL properties in common with 
the theory of HLR, including a finite compressibility, yet 
differs in that there is no CS term for the gauge field, 
while the physical density is dipolar or (using an equa- 
tion of motion) is — pV A a. 

In the rest of the paper, we follow a different argument 
from that just presented. We give a detailed microscopic 
derivation, in which the relationship p(r) = — pSJ A a ap- 
pears only at the end; thus we do not rely on the Berry 
phase argument. The starting point is an approach of 
Pasquier and Haldane, described in Sec. Ill below. In 
this approach, which works for cf> — I only, that is bosons 
at v = 1, each fermion is described by two coordinates, 
which we term "left" and "right" , but the available states 
are those of a particle in zero magnetic field, because the 
wavefunctions are complex analytic in the left and anti- 
analytic in the right coordinates. The left coordinate is 
that of the underlying particle contained in the fermion, 
while the right coordinate rep resen ts an attached vor- 
as in the pictures in Sec. II B. The system must 

'' " in the right 



tex, 

obey a constraint of fixed density p ±l = p 
coordinates. Since the separation of the left from the 
right coordinate is Ak when the fermion is in a plane 
wave state of wavevector k, the physical density is dipo- 
lar. In order to maintain the constraint, the longitudinal 
part of the current j R of the vortices (right coordinates) 
must vanish, as argued above. In Sec. we consider a 
conserving approximation for observable response func- 
tions. We show that the constraints are satisfied in this 
method. We calculate the density-density response, its 
spectral density, the longitudinal conductivity, the scat- 
tering of the fermions by a potential, and the current- 
density operator. From the results we deduce that the 
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system can be described in terms of the strongly-coupled 
gauge field mentioned above. The gauge invariance is a 
manifestation of the constraint. The gauge fields obey 
the CS equations, even though there is no CS term in 
the action. In Sec. we indicate the form we expect 
for the exact results to all orders in the interactions, and 
give arguments why these are correct. We conjecture 
that a certain sum rule for the spectral density is exact. 
While at present this approach works for <fi = 1, that 
is for bosons at v = 1, we expect that the conclusions 
are more general, as the results and arguments of the 
previous subsections and the beginning of this one are. 

III. PASQUIER-HALDANE APPROACH 

FOR cj> = 1 

In this section we review (with a few variations of our 
own) the method of PH J52| which works only for <f> — 1 , 
though the filling factor does not necessarily have to be 
one. A similar method works for fermions with one vor- 
tex attached, mapping them to composite bosons. Since 
the formalism has not appeared elsewhere in the form in 
which we will use it, it will be presented in self-contained 
fashion. 

We begin abstractly, labelling arbitrary single-particle 
states with indices. Hopefully the later development in 
coordinate space, though less general, will seem less ab- 
stract and give more physical insight, and show clearly 
the connection with composites particles and the LLL. 
We take fermion operators which are matrices with two 
indices, c mn and c\ lm , with canonical anticommutation 
relations 

{c mn ,C n / m /} — &mm'&nn' (3T) 

(and others vanish) where m, m! , n, n' run from 1 to N 
(this case of square matrices is convenient for the v = 1 
boson problem, while rectangular matrices would be used 
for v ^ 1). The anticommutation relations arc invariant 
under independent unitary transformations on the left 
and right indices, under which 

c i-> U L cU R , 

C t i ^ Uyul, (3.2) 

where Ul, Ur are unitary N x N matrices. These trans- 
formations are generated by the operators 

Pnn' ~ ^ ] c nm c mn' , (3-3) 
m 

Prnm' ~ ^ ] C nm' Cmn - (3-4) 
n 

The right generators p R generate the group \J(N)r of 
unitary matrices. These are used to specify a set of iV 2 
constraints on the system, 



0C' - <W)|*phys> = (3.5) 

which defines a subspace of states that will be iden- 
tified with the physical Hilbert space. By taking the 
trace, we see that these imply that the U(l) generator or 
fermion number operator (which is common to \J(N)r 
and U{N) L ) 

N = J2c f nm c mn , (3.6) 

mn 

must have eigenvalue equal to N. Thus in the allowed 
subspace, N is both the range of the indices, and the 
number of fermions. The remaining right generators gen- 
erate S\J(N)r, and physical states must be singlets un- 
der the action of this group. The other group, SXJ(N)l, 
is not used for constraints, and will be broken by the 
Hamiltonian to a subgroup that represents translations 
and/or rotations on the two-dimensional manifold (say, 
the sphere, torus, or infinite plane) on which the physical 
particles move. At the same time, the generators p^ m / 
will represent the physical density on this manifold. 

The physical states that satisfy the constraints can be 
written as linear combinations of 

i^ r mi...m N \ _ 7U...7WV t c t ... c t |Q\ 

n\ ,. . . ,njy 

(3.7) 

where |0) is the vacuum containing no fermions. These 
states contain N fermions and are clearly singlets under 
SXJ(N)r since they are antisymmetric in the n (right) 
indices. On the other hand, the anticommutation of the 
cT's implies that they are symmetric in the remaining m 
(left) indices. Thus these states can be viewed as basis 
states for a system of N bosons, each of which can be in 
any one of N single-particle states. Such a boson system 
could be described by basis states 

aL-'-aljO), (3.8) 

where [a m ,o^ m ,] — S mm > and others vanish. Each such 
state is obtained in this way, which proves that the 
fermion system of c's with the constraints is equivalent 
to the unconstrained boson system. If we define a filling 
factor as the particle number divided by the number of 
available orbitals, as N — > oo, then in our case we clearly 
have bosons at filling factor v = 1. 

We note that in the larger Hilbert space without the 
constraints, which is just the Fock space of the c's, each 
fermion can be in any of iV 2 states, so there are 




linearly-independent states for N fermions. The states 
satisfying the constraints form the Fock space of the 
bosons a, which contains only 
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2N- 1 
N 



(3.10) 



linearly-independent states. 

The left indices m can run over any range, and this can 
be used to represent any filling factor v. The constrained 
system can also be set up using canonical commutation 
relations for the c's, and a similar argument then shows 
that the physical states represent fermions, (e.g. elec- 
trons) at v < 1. 

So far we have only a way of representing bosons by 
fermions (or vice versa), and the technique is reminis- 
cent of the methods used for quantum spin systems, in 
the case of a single quantum spin (see e.g. 53 1). If we now 
view m, n as indices for lowest Landau level states, say 
on a sphere where there are Nj,+ 1 such states for flux 
quanta through the sphere Eq| , then for the case where 
both indices range from 1 to N, we have N = N<f, + 1, 
and the filling factor v agrees with that defined as N/N^ 
as N — > oo. We can introduce coordinate space wave- 
functions for the left index m, which are just those of the 
physical bosons. We do the same for the right indices 
n, except that they are complex conjugated so that the 
field strength (or the charge) is effectively reversed. Us- 
ing orthonormal single-particle LLL basis states u m (z), 
we write in analogy with the usual field operators 



c(z,w) = };U m (z)u n (w)cm n , 

mn 



(3-11) 



which are adjoints of each other. Note that we use z's 
for "left" indices, corresponding to m's (which however 
appear on the right in c') and w's for "right" indices, 
corresponding to n's. The appearance of two coordi- 
nates on c and means that they behave like opera- 
tors on the LLL single-particle Hilbert space, just like 
the matrix structure they had in index notation. A for- 
malism for handling such operators as integral kernels 
is given in Appendix |A[ For the sphere, we can write 
u m (z) oc z m , for m — 0, . . . , = N — 1, and the factor 
(1 + |z| 2 /4_R 2 )^ (JV * +2)/2 must be attached before integra- 
tion. Following this convention we will write only the 
polynomial part in the following wavefunctions. 
In the z, w variables, the densities become 



p R (w,w') = / d 2 zc^(w, z)c(z, w/), 



o L (z,z') 



d 2 w c'(w, z')c(z, w) 



(3.12) 
(3.13) 



Matrix multiplication has been replaced by integration, 
so that all operators in the single-particle Hilbert space 
of LLL functions of z and W become integral kernels (see 
Appendix [A]). One can see that p L (z,z) is the LLL- 
projected density operator denoted p by Girvin, Mac- 
Donald and Platzman (GMP) [p0[, and p R is analogous. 



Passing to the thermodynamic limit at fixed field 
strength and density = p, the radius of the sphere goes to 
infinity, the system becomes flat locally, and we may use 
Fourier transforms. The version of the Fourier transform 
required is defined in Appendix |A]. To avoid discussion 
of global issues, which would distinguish this thermody- 
namic limit from that of a torus, we will view the use 
of Fourier transforms as a technique for handling local 
calculations, in which we could include damping factors 
which tend to unity at the end. Alternatively, every cal- 
culation could, with only a little extra difficulty, be done 
in coordinate space. A third alternative would be to use 
the analog of the Fourier transform, involving spherical 
harmonics, on the finite size sphere. This is more tedious. 
Introducing the Fourier transform in the plane, then, we 
notice that the pair of coordinates z, W for each particle 
or field operator c is replaced by a single ordinary two- 
dimensional wavevector k. This makes sense because, by 
choosing equal and opposite field strengths for the basis 
functions in these coordinates, the particles effectively 
see zero magnetic field, for our filling factor v = l/(f> = 1. 
Note that, because the functions are analytic in z, w (the 
LLL restriction) , we do not have effectively four real vari- 
ables per particle, as we would if the basis states had not 
been restricted to the LLL. The transformation of the 
matrix c(z,w) into a plane wave operator is similar to 
that for the density operator, say p L , which can clearly 
be traded for its Fourier components (see e.g. GMP). 

In terms of Ck, c^, which are defined in Appendix |A|, 
and which satisfy 



{c k ,cj c ,} = (2 7 r) 2 5(k-k'), 



we have 

p r (<i) 



d 2 k 

(27T) 2 
d 2 k 



a -iikAq t 



k-iq C k+iq> 



and we can show that 
[/Aq),p L (q')] =0, 



,iikAq t 



k l q C k+lq, 



(3.14) 

(3.15) 
(3.16) 

(3.17) 



[p fl (q),p fl (q')] =-2isin-qAqV(q + q'), (3.18) 
[p L (q),P L (q')] =2isiniqAqV(q + q')- (3.19) 

The L ie algebra commutation relations defined by Eq. 
(|3.19|) appeared in GMP and in Ref. @, and the alge- 
bra so-defined has become known as (the defining 
relations are often given in a different basis of the Lie 
algebra, essentially the expansion of our p L (z,z') in an- 
gular momentum eigenstates z m , z' m ). In the notation 



of GMP, our p L (q) 



p(q). (The following alge- 



braic comments will not be used in the following.) From 
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our point of view, Woo is just a certain limit of SU(iV") as 
N — > oo. It is also helpful to note that if the 2 sin ^qAq' is 
replaced by qAq' in the commutation relations (for exam- 
ple, because q, q' or the magnetic length are small), then 
the resulting algebra is that of "area-preserving diffco- 
morphisms" , or equivalently (for the corresponding Pois- 
son bracket relations) Fourier components of functions 
on classical phase space. Woo can then be viewed as a 
quantum deformation of the latter, thus as "diffeomor- 
phisms of the quantum analogue of phase space" , a fairly 
familiar view of the LLL. The connection of W>o with the 
quantum Hall effect has often been remarked |55| ]. Our 
interest here is in the isomorphic algebra generated by 
the p R, s, which are the constraints of our problem. 
The constraints become 



(p R (q) - p(2n) 2 S(q))\^ phys ) = 0. 



(3.20) 



Thus states can be built up in the "big" Hilbert space as 
combinations of 



n 4io) 



(3.21) 



(where the product is indexed by k's in a set of N 
wavevectors kj), and then projected to satisfy the con- 
straints. The effect of projection can be more easily ap- 
preciated in terms of wavefunctions in coordinate space, 
by returning to the finite size system. 

In coordinate space, the constraints require that the w 
dependence of wavefunctions be that of a full LLL, 

*phys(^l) Wl, ■ ■ - ,ZN,W N ) = f(zt, ...,Zn) J\(Wi 



(3.22) 



because the LJ factor in the W's is the unique totally- 
antisymmetric function annihilated by the p s, since the 
full LLL has no density fluctuations. Hence / is a sym- 
metric polynomial in the Zi's, as appropriate for bosons. 
Projection of the wavefunction of any state in the "big" 
Hilbert space to this physical subspace, where states can 
be characterized just by /, is accomplished by multiply- 
ing by Jl^i — w j) an d integrating over the w^s with 
the appropriate measure, leaving a symmetric function 
/ in the z^s (possibly zero). If as a family of exam- 
ples we take the states ( |3.2l[ ), or their analogues on the 
sphere, in first quantization they become Slater determi- 
nants det(lL i M i (zj, Wj)), where the Ylm(z, w) are spher- 
ical harmonics projected to the LLL, which correspond 
to the plane waves Tk in the plane, defined in Appendix 
[a|. Then the projection gives 

/ = / II d2u ' k Y\( W i ~ W j) det Y L,hh ( z j,Wj) 
k i,j 

= 7>lll det Y LiMi (%) Y[(z t - Zj ), (3.23) 



that is, the projection to the LLL of ordinary spherical 
harmonics in a Slater determinant times the LJ factor. 
The se are just the trial wavefunctions described in Sec. 
|Tb|. Thus the formalism not only describes bosons at 
v = 1, but the fermions are closely related to those in 
the "physical" approach, where the amplitude of the LJ 
factor is automatically included in the trial wavefunc- 
tions. Contrast this with the CS approach, where the 
trial wavefunctions satisfying the CS constraint of one 
flux attached to each particle consist of the Slater de- 
terminant times only the phase of the LJ factor, and no 
LLL projection. Note also that while the projection into 
a strictly smaller subspace implies that states described 
by distinct sets of k^ before projection may not be orthog- 
onal after projection, they do not usually vanish, except 
in some exceptional cases noted in Ref. [Q . 

Since the right coordinates w of the fermions become, 
in the trial wavefunctions after projection, the locations 
of the vortices, it seems natural to refer to them as 
such even before projection. Thus we can say that each 
fermion consists of a particle (boson) at the left coordi- 
nate z, and a vortex at the right coordinate w, and so 
as a whole is effectively neutral. The constraints demand 
that the density p R of vortex coordinates is fixed, as an 
operator statement. This seems natural if the vortices 
are thought of as forming a two-dimensional plasma (in 
view of the LJ factor and Laughlin's plasma mapping 
]l4||), since the plasma is in a screening phase and sup- 
presses long- wavelength density fluctuations; indeed, in 
this case of v = 1, there are no fluctuations in the LLL 
density at all in the Laughlin state (the full LLL or Van- 
demonde determinant). In retrospect, this condition on 
the vortices seems to be the main effect that was left out 
in Refs. 

Now we finally specify the Hamiltonian appropriate to 
bosons in the LLL at v = 1. In terms of the boson oper- 
ators a introduced earlier, we have, assuming a potential 
interaction between the bosons, 



2 ^ 

mi ,...,7714 



v, 



at at 



mini2;m3m4 u m 1 Ll mo ^ m 4 u m3 i 



. (I , . 



(3.24) 



where the matrix elements of the interaction in the LLL 
are |5q] 



V, 



mim2;ni3m4 



d 2 ri d 2 r 2 u mi (zi)u m2 (z 2 ) 
xV(ri - r 2 )u m3 (z 1 )u mi {z 2 ). (3.25) 



The corresponding operator in the large Hilbert space, 
where it commutes with the constraints p R , and so 



projects to H in Eq. (3.24), is 



1 



H — ^ y ] Vmim 2 \m a m i cl limi cl l2m2 C m4n2 C 7n3ni . 
mi ,...,m4 



(3.26) 
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Then using the definition of c{z, w) we obtain 

H =\J d 2 r 1 d 2 r 2 V{r l - r 2 ) : p L (z 1 ,z 1 )p L (z 2 ,z 2 ) :, 

(3.27) 

where the normal ordering is with respect to the vacuum 
of the c's, |0). Thus this is simply a potential interac- 
tion written in terms of the LLL-projected density p L . 
In Fourier space this becomes 

h = \J : p L (q)p L (-q) =. (3-28) 

where V(q) — e^^ q ^V(q) absorbs a factor left from the 
definition of the Fourier transform of p L , and V(q) is the 
usual Fourier transform 



V(q)= J d 2 re~ tcl r V(r) 



(3.29) 



The interaction Hamiltonian breaks the symmetry group 
from SU(N)l (in the absence of interacton) to SU(2) (for 
the sphere) or to magnetic translations and rotations in 
the case of the plane. It still commutes with the "con- 
straint operators" 



G(q) = p R (q) 



p(2n) 2 6(q). 



(3.30) 



Our expression for the Hamiltonian differs somewhat 
from that in the paper of PH. They work on the torus, 
which is a relatively unimportant difference, and write 
the Hamiltonian using the constraints to make the ansatz 
explained in Sec. II C , which results in a one-body term 
that gives the fermions an effective kinetic energy coming 
from the interaction. In our approach we do not wish to 
make such a substitution since the commutator of H with 
the G(q) would not vanish identically, but only on using 
the conditions G(q) = 0. The reason for our insistence 
on retaining [H, G(q)] = will be discussed in the next 
section. Of course, if everything is done correctly, the 
results should be the same, in the end, since the starting 
point is the same. 



IV. HARTREE-FOCK AND CONSERVING 
APPROXIMATIONS 

In this section, which is the central one of the paper, 
we develop an approximate solution for our system that 
descibes the FL state. We begin in Subsec. IV A with 
the Hartree-Fock approximation, which yields a disper- 
sion relation for the fermions. Then in Subsec. 
we explain how the constraints can be included 



IV B 



We 

choose a gauge such that for nonzero frequencies, they 
must be satisfied without any assistance from integration 
over auxiliary fields that impose them explicitly. This is 



achieved in Subsec. IV C by use of conserving approx- 
imations, a familiar method of many-body and quan- 
tum field theory. In the present case, such an approx- 
imation consistent with HF is the generalized or time- 
dependent HF approximation, which sums ring and lad- 
der diagrams. We show explicitly that the constraints are 
obeyed in our approximation. In Subsec. IV D| we inves- 
tigate the asymptotics of the ladder series that appears 
in Subsec. IV C, for use in the following calculations. In 



Subsec. IV E we apply the approach to the physical re- 



sponse functions, beginning with the density-density re- 
sponse. We show that the system is compressible and 
that the longitudinal conductivity relevant for the sur- 
face acoustic wave experiments, which is a certain limit 
of this response, is given by exactly the same expression 
as in HLR. We also exhibit a sum-rule-like relation for 
the high frequency response, or for the first moment of 
the spectral density, which we will later argue is exact. 
We consider the scattering of a fermion by a scalar po- 
tential perturbation, and interpret the result in terms of 
a vector potential related to the density by the CS rela- 
tion discussed in Sec. 0. We calculate the longitudinal 
conductivity due to impurity scattering. Finally, we con- 
sider the physical current density, which we relate to the 
stress or momentum flux tensor of the fermions, and so 
recover the other CS relation. 



A. Hartree-Fock approximation 

In this subsection, we use the HF approximation, 
which is quick and is the simplest one that gives an effec- 
tive kinetic energy and is consistent with a stable Fermi 
sea as the ground state. The treatment of the constraints 
will be extensively discussed in the next subsection, and 
the formalization of the exchange part of the self energy 
as the saddle point approximation to a functional inte- 
gral, valid in some sense in a large-M limit (in a gener- 
alization of the model to M component fermions), is left 
to Appendix |b]. 

The problem for tf> = 1 us ing t he PH approach is de- 
scribed by the Hamiltonian (3.28) which can be written 



d 2 hd 2 k 2 d*q ~ UklAQ _ 
2 J (2tt) 6 W 



"ki- 



, c k 2 -iq c ki + iq> 



(4.1) 



subject to the constraints G(q) = p R (q) — p(2ir) 2 S(q) = 
0, that is N = N, and 



d 2 k 
(2^) 2 



ikAq t 
°k- 



-, C k4 



= 



(4.2) 



for q 7^ 0. Notice that when the phase factor containing 
k A q is expanded in a Taylor series, to 0(q 2 ) it takes 
the same form as the constraint found by SM and Lee 
p4|,E5||, as mentioned in Sec. |IIC. 
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The HF approximation for a translationally-invariant 
system takes the energy eigenstates to be Slater deter- 
minants of plane waves, that is plane-wave-occupation- 
number eigenstates in the second-quantized formalism 
and the energy of such a state to be the expectation value 
of H . As is well-known, for the excitation spectrum, this 
is equivalent to replacing H by an effective one-body 
Hamiltonian with an effective energy £k for each plane 
wave state k, where £k depends self-consistently on the 
occupation numbers n^. In the present case, we must 
also include the constraints by the use of Lagrange mul- 
tipliers A q and minimise 



H-fiN- 



d 2 q 
(2tt) 2 



A q G(-q) 



(4.3) 



with respect to A q to find the ground state. When al- 
most all particles are in the Fermi sea, the A q are zero 
by translational symmetry, except at q = where Ao ab- 
sorbs the chemical potential fi, consistent with the fact 
that the constraints fix the particle number and hence 
we are actually in the canonical, not grand canonical, 
ensemble. Consequently one has A q = (27r) 2 A<5(q), and 
A + /i is determined by the condition on the total par- 
ticle number. One arrives therefore at the total energy 
expectation value, 



kk' 



(4.4) 



(in which we have used the conventional notation for a 
finite system in a square box of side L, with discrete k 
values, and 7ik are the expectation values of the occupa- 
tion numbers for the corresponding states), where 



/kk' = ^(0)-y(k-k'). 



(4.5) 



The function /kk' plays the role of the Landau inter- 
action function when k and k' are restricted to the 
Fermi surface. The effective single-particle Hamiltonian 
K = H-(fi+ X)N is 



k 



£kC k C k , 



where £k = £k — A 4 — ^ an d 



£k = V(0) 



d 2 k> 
(2ir) 2 



d 2 V 



k') 



l k'j 



(4.6) 



(4.7) 



in which the first term is the direct or Hartree term, equal 
to V(0)p, and the second is the exchange or Fock term, 
which is responsible for the k-dependence of £k- Also, in 
the ground state at zero temperature, n£ = d(kp — k), 
kp = v2 in our units, and /x+A is chosen so that £fc F = 0. 
Notice that the phase factors in the Hamiltonian H have 
turned out to be unity in the HF expressions, which are 
identical to those of the usual Fermi gas, except that the 



bare kinetic energy is zero, and that V"(q) replaces V(q) 
for reasons connected with the LLL. This formula for £k 
differs from that of other authors, discussed in Sec. [II Cj , 
in that it depends explicitly on the occupation numbers 
of the other k states, and does not reduce to the self inter- 
action of a dipole even for small q = k k' in the integral 
in the exchange term. Our £k gets its k dependence from 
the exchange effect, while the interaction of the parti- 
cle with the correlation hole that surrounds it (due to 
the vortices) is a "Hartree-likc" term (and not simple 
Hartree) (see (l(J where exchange effects were explicitly 
neglected). Thus the exchange effect found here in the 
simplest approximation seems to be complementary to 
the interaction with the correlation hole, and probably 
both terms would be present in a better approximation. 
As for the dipolar form of density, we will see that the 
density does take on this form, and this could be in- 
cluded in the exchange self energy, but this would ne- 
cessitate a complicated self-consistent calculation which 
could not be done analytically. In any case, the dipolar 
effect changes the form of the interaction at small q, while 
intermediate q values are important in the exchange self 
energy. Thus the expression here is a convenient starting 
point, and not badly wrong physically, at least in some 
cases, as we will see shortly. 

The zero-temperature HF dispersion relation can be 
studied in detail. Apparently, no difficulties are caused 
by the absence of a bare £k term. For any repulsive 
interaction V(q) = e~2' q ' V(q) > 0, £k increases mono- 
tonically with |k| for all k. At |k| = kp, 



kp 
m* 



d\k\ 



d 2 k dV 

i k '| <fcF (2^^kl 

kp 
2^ 



(k'-k) 



^lF(k'-k)cos0 kk , 



(4.8) 



(note that #kk' parametrises the angle between k' and 
k which are both on the Fermi surface). For a 8- 
function (short-range) potential, V(q) = V(0), l/m* is 
positive and finite. Thus the system is stable against 
single-particle excitations. For a Coulomb interaction, 
V(q) = 27re 2 /|q|, there is a logarithmic singularity at 
|k| = k F : 



d£k 
d|k| 



lnlfe 



(4.9) 



This is very similar to that for the Coulomb interaction 
in the three dimensional electron gas at zero magnetic 
field treated in HF approximation. In that case, the di- 
vergence is unphysical and is removed by replacing the 
bare Coulomb interaction in the exchange term by the 
screened one, which leaves a finite effective mass and heat 
capacity Cy ~ "fT ~ m*kpT. This conclusion of course 
depends on the presence of screening due to the nonzero 
compressiblity of the electron gas. In the present prob- 
lem, the existence of such a compressiblity is one of the 
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points we wish to study, so we must return to this later. 
Note, however, that replacing the unscreened interaction 
by the dipolar interaction also cuts off the divergence 
in the present problem. As mentioned already, this will 
also be left for later discussion. For the time being, we 
may consider an interaction of shorter-range (decaying 
as a faster power) than Coulomb and the effective mass 
is then finite within HF. 

The question may be raised of whether a charge- 
density wave (CDW) instability could take place due 
to the absence of a bare kinetic energy. However, 
the constraints p R (q) = 0, though not the same as 
J d 2 k i q c k+ i q = 0, may have a similar effect in main- 
taining the uniform density of the fluid within HF (a 
CDW in the underlying particles cannot be ruled out at 
some filling factors, especially but may not be de- 

scribable within HF for the fermions). Another possible 
instability is to pairing as in BCS theory. This has been 
argued by PH p7| ], who found numerically that bosons 
at v = 1 tend to form a ground state with high over- 
lap with the Pfaffian state, a paired state which is pre- 
sumably incompressible. However, for some interactions, 
such pairing may either not occur, or be very weak so that 
it occurs only at very low energies, and then the present 
results for the "normal" Fermi-liquid-like state will still 
apply at higher energies, temperatures, or wavevectors. 
For the state of electrons at v = 1/2, experiment and nu- 
merical results both indicate that pairing must be either 
extremely weak or absent, so there would seem to be a 
regime to which the theory would apply, assuming that 
it can be extended to (f> > 1. We return to the issue of 
pairing in Sec. 0. 

B. Constraints 

In this subsection we begin a fuller and more system- 
atic analysis which begins from the HF approximation 
but entails a careful study of the role of the constraints. 
In the present subsection, we explain a functional inte- 
gral method for handling the constraints exactly. Ap- 
proximation methods are discussed beginning in the fol- 
lowing subsection, where the starting point is once again 
HF. The present subsection could be skipped on a first 
reading, but does explain why many statements later in 
the paper are restricted to nonzero frequencies. 

The constraint operators G(q) obey 

[G(q), G(q')] = -G(q + q')2i sin lq A q', (4.10) 
[#,G(q)]=0. (4.11) 

These relations have the property that if all G(q) are 
replaced by zero throughout, as stipulated by the con- 
straint, then they are still true. Constraints with this 



property are termed first-class, while others are termed 
second class pSTj] . Second-class constraints lead to modi- 
fied commutation relations given by "Dirac brackets" in 
the constrained subspace, and are generally more awk- 
ward to handle. An example is the constraint of being 
in the LLL, applied to one or more charged particles in a 
magnetic field, which when imposed in the obvious way 
is second class, and consequently the coordinates x, y 
of the particle(s) end up not commuting when projected 
into the LLL. By contrast, systems with only first-class 
constraints can be viewed as gauge theories and there 
are very well-developed methods by which they can be 
handled f57fl . The advantage of the PH approach is that, 
while the fields are in the LLL from the beginning, the 
only constraints involved are first class. 

The importance of the first-class property of the con- 
straints is that G(q) form a Lie algebra, SU(iV) or Woo, 
and are constants of the motion, dG(q)/dt = for all q. 
Thus, before considering them as constraints, the G(q) 
can be viewed as generators of a symmetry algebra of 
the Hamiltonian. As constants of the motion, the condi- 
tions G(q) = 0, if imposed at the initial time, would hold 
for all other times. Our procedure, which is a version of 
the Faddev-Popov functional integral method, will differ 
somewhat from this, however. To find thermodynamic 
properties and correlation functions, we begin with the 
partition function, 

Z = Tr G=0 e~^ H -^\ (4.12) 

where the trace is restricted to states satisfying the con- 
straints. This can be written formally as 

Z = TYe-^ H -^ ) 5 Gfil (4.13) 

where the trace is taken in the Hilbert space, the Fock 
space of the fermions c, with no restriction on the fermion 
number N . (The /iN term is included to make this look 
conventional, even though the constraints fix N = N, so 
the constrained ensemble is canonical, not grand canon- 
ical.) The 8- function, which imposes all the constraints, 
can be given a Fourier representation which essentially, 
for a nonabelian group, means integration over the group 
manifold. Here we return to the U(A) notation that we 
had for finite N: 

§G,o = J [U-'dU] U fj g e^- N \ (4.14) 

where the first integration is over SU(A) with the invari- 
ant (Haar) normalized measure [U~ 1 dU], and the sec- 
ond is over U(l) and imposes N — N. We can write 
U = e ^2 a l f )x " G " ( w here a = 1, . . . , N 2 — 1 runs over 
a basis of the SU(A) Lie algebra) and convert the unre- 
stricted Tr to a functional integral in the standard way 
to obtain 
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z = 



V[c,^][U- l dU]^^- exp 

Z7T 



dr < Tr — c - 



i ^ ' A a G a 



iXo(N-N) 



(4.15) 



where i?, TV, and G a are given by the standard forms in terms of the Grassman variables c mn (r), cjj m (r), and the trace 
in the exponent is on the U(JV) indices. The commutation properties (4.11) were used in obtaining this expression. 
The A a 's and Ao = 9/(3 now play the role of time-independent scalar potentials in the sense of gauge theory. The 
functional integral results from gauge fixing a manifestly gauge-invariant version, 



T>[c,J]T>[<j)] exp 



r/r 



. d 



)c f c- 



H-fiN 



(4.16) 



in which cf> stands for all the A's in N x N matrix form, 
is r-dependent, and is functionally integrated over the 
U(iV) Lie algebra. Under a U(JV) gauge transformation 
U,(j)^> U~ x <j)U + U^dU/dr. This reduces to the previ- 



ous integral (4.15) by imposing the condition d(f>/dr = 
inside the functional integral (we are neglecting Faddeev- 
Popov determinants). This condition is not the same as 
= (which is often used instead), which cannot be 
reached by a gauge transformation from an arbitrary </>, 
since gauge transformations must be periodic in r with 
period (3. Thus J dr <f> cannot be gauged away to zero. 

The holonomy Pel dT ^ (P denotes that the integral is 
path ordered), which is an element of the group U(JV), re- 
mains. This holonomy is the combination U e lS of the ear- 
lier integration variables. Under a r- independent gauge 
transformation it is not invariant: 



Pe 



J dr<t> 



V- x Pe 



J dT<f> 



u. 



(4.17) 



and so only the set of eigenvalues of this matrix is gauge 
invariant. (Note that there are gauge transformations 
that permute the eigenvalues.) The integral in eq. ( 4.15 ) 
is over the holonomy, but can be further gauge-fixed to 
leave integration over the eigenvalues only: 



J [u^du]^!^^^ 



Nl 



xe 



a=l 1 ' 7 <5 



(4.18) 



with the measure well-known in, for example, random 
matrix theory (which here has no connection with the 
similar-looking LJ factors!). 



The reduction of the constraint integrals to only zero- 
frequency fields shows that at low temperatures, the inte- 
gration over these fields is relatively unimportant, since 
zero frequency is of zero measure in integrals over fre- 
quency that appear in a diagrammatic treatment, as will 
be used in the following. The non-zero frequency part of 
the constraints G(q, ui) = will have to come out auto- 
matically without help from an integration over a field 
that enforce s it d irectly (as in the totally gauge unfixed 
version eq. ( 4.16| )). It will be demonstrated that this oc- 
curs in the next subsection. 

Finally we note that when developing the HF approxi- 



mation as in Sec. IV A (or when taking the saddle point of 
the functional integral as in Appendix |b|), the Lagrange 
multiplier A is the saddle point value of iAo, so the sad- 
dle point value of Ao is imaginary. This phenomenon is 
common in such treatments. 



C. Conserving approximations 

In this subsection we ret urn to the approximate treat- 
ment begun in Sec. EVA, consider response functions, 
and address the question of whether the constraints are 
satisfied. The central issue is the use of a so-called con- 
serving approximation, that is an approximation that 
satisfies the relevant Ward identities, which express the 
symmetry under V(N) or generated by the con- 
straint operators G(q). 

The appropriate conserving approximation to use for, 
say, the density-density response in a normal Fermi liq- 
uid depends on the approximation used for the one- 
particle properties, that is, the conserving property in- 
volves consistency of approximations for different prop- 
erties. It is well-known that the random-phase approx- 
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imation (RPA) corresponds in this sense to the Hartree 
approximation, and perhaps less well-known that gener- 
alized RPA, also called time-dependent HF, corresponds 
to the HF approximation (for discussion of conserving 
approximations, see e.g. [ |58|j59| ]; for the generalized HF 
approximation in a FL, see PN, Ch. 5). These are some- 
times stated in terms of $-derivability, that is approxi- 
mations that can be derived by making an approximation 
once and for all for the free energy <I> (or for the thermo- 
dynamic potential) in the presence of source fields that 
couple to the onservables of interest (such as the den- 
sity), and then obtaining response functions in the same 
approximation by taking functional derivatives with re- 
spect to the sources, guaranteeing the same sort of con- 
sistency. 

The importance of the conserving approximation dep- 
nds on the nature of the problem. In the example of a 
normal Fermi liquid, the basic symmetry is conservation 
of total particle number, which is not broken by Hartree 
or HF. The conserving approximation is then needed to 
ensure that the Fermi liquid relations are satisfied, pro- 
viding detailed relations among physical quantities. By 
contrast, in a BCS superconductor, the simplest approx- 
imation (which can be viewed as an extension of HF) vi- 
olates conservation of particle number, and the conserv- 
ing approximation |58[] not only restores gauge invariance 
(number conservation) but also leads to the prediction of 
a collective mode, the Anderson-Bogoliubov mode (which 
is the Goldstone mode connected with the spontaneous 
symmetry breaking in the case of short-range interac- 
tions). Thus the use of a correct approximation has much 
greater physical consequences in the latter case. 

Turning now to the present problem, the HF approxi- 



mation of Sec. IV A does not break conservation of total 



particle number N. However, the symmetry generators 
G(q) = p R (q) for q ^ are not conserved by the HF 
approximation as it stands. The easiest way to see this 
is that G(q) does not annihilate the HF ground state, 
which is just the Fermi sea \FS). Thus this state does 
not satisfy the constraints G(q)\FS) = for q ^ 0. It i s 
also clear that the HF effective Hamiltonian, Eq, (4.6), 
does not commute with these G(q). The solution to this 
problem will have to use the conserving approximation 
appropriate to our HF approximation. Since there is a 
conserved quantity for all q, the results will be even more 
striking than in cases such as the BCS theory where only 
a global symmetry was broken. We note that the Fermi 
sea can be made invariant by projecting to the invar- 
ian subspace as in Eq. (3.23). However, such a projec- 
tion necessitates that further work be numerical. Ana- 
lytical work, and thus conceptual understanding, can be 
achieved only by persevering with the gauge theory ap- 
proach. Rather than give up the Fermi sea trial state and 
the HF energies and searching for some other, invariant, 
starting point, we keep it and take care of the constraints 
by the following conserving approximation. 



The conserving approximation will be illustrated here 
by the calculation of the p R -p R , p R -p L and p L -p L 
imaginary-time response functions (more precisely, the 
generalized susceptibilities), defined in Fourier space by 



X lJ (q,w„)(27r) 2 (5(q + q')/3(S w „+^„,,o = 

(p l (q^n)p 3 (q'^n0>, (4-19) 

in which i, j can be R or L, tu n are the usual Matsubara 
frequencies, and it is implicit that the connected part of 
the function is taken, thus dropping a ^-function term 
containing (/9 l )'s. The conserving approximation that 
corresponds to HF takes the form of the sum of all ring 
and ladder diagrams. The Green's function lines in the 
diagrams are the HF Green's functions 



(4.20) 



The usual Dyson-equation argument leads to formulas 
in terms of the one- inter action irreducible susceptibilies, 
as discussed in Sec. HA , defined as those diagrams that 
do not become disconnected when one interaction line 
is cut (note that we disregard the Hartree self-energy 
diagrams that are implicitly included in out HF Green's 
functions, which means we are treating the diagrams here 
as skeleton diagrams; such terms would be absent anyway 
for a long-range interaction due to the neutralising back- 
ground). These formulas, which are completely general, 
are (all %'s have the same arguments q, ui n ) 



Xll 



Xrl 



irr 

ALL 



1 + V(q)xf L 



Xrl 



1 + V(q)xf L 



XRR = Xrr ~ Xrl 



1 + V(q) X f L 



irr 

Xlr- 



(4.21) 
(4.22) 
(4.23) 



Note also that XLi?(q,w„) = Xrl(~<1, -w n ). The con- 
serving approximation is now the statement that the var- 
ious x' rr are to be calculated (for w ^ 0) as the sum 
of the ladder diagrams, with the HF Green's functions. 
Since p L is the physical density, x'll ^ s the one of most 
physical interest for long-range V^(q), such as Coulomb 
interactions. 

We begin with Xrri so as to show that at w ^ 0, 
the fluctuations in the constraints G(q) vanish in our 
approximation. The Feynman rule for the interaction 
can be read off in the standard way (5(| ; it includes the 
wavevector-dependent phase factor as well as V(q). Also, 
there is a phase factor in the p R vertices, as in eq. ( 3.15 ) . 
Note that those in the interaction arise from the phase 
factors in the physical density p L , Eq. (3.16). In the lad- 
der diagrams for x'rr the structure of the momenta is 
such that all the phase factors cancel, as the industrious 
reader will verify. Note that this is an exact statement, 
and not only valid at small wavevectors, whether internal 
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or external, so the exponential defining the phase factor 
was not expanded in a Taylor series. Consequently, for 
the ladder diagrams for X^r only, the ladder series is 
identical to the same approximation to the irreducible 
susceptibility in the usual density 



p(q) 



d 2 k 

(27T) 



2 L k- 



(4.24) 



in a model with Hamiltonian 

1 f d 2 k, d 2 k 2 d?q - . 



(4.25) 



with no kinetic-energy term. This could be phrased by 
saying that there is the ordinary, Galilean-invariant ki- 
netic energy term with zero magnetic field, but the mass 
mo is infinite. We call this latter model the zero field, 
infinite mass (ZFIM) model. Note that the HF approx- 
imations in the two models also coincide, because the 
phase factors disappeared there also. In the ZFIM model, 
[p(q),H] = for all q, so the model possesses a gauge 



symmetry, whether or not we wish to impose a constraint 
p = constant. In fact, if such a constraint were imposed 
in this model, there would be no states that satisfied it at 
all. The reason (in classical language) is that in a contin- 
uum model, any configuration of point particles clearly 
has nonconstant density. In a similar model on a lattice, 
solutions to the constraint exist only if the value of the 
particle number required by the constraint at each site is 
an integer, since these are the eigenvalues of the number 
operator for each site. This cannot be satisfied if we take 
the continuum limit (zero lattice spacing) at fixed aver- 
age density. In our system representing the LLL, which 
is in the continuum, many solutions to the constraint do 
exist, provided we choose (similarly to the lattice ZFIM 
model) the constrained value of the total number to be 
the same as the range of the right indices n, as we have 
done. Therefore, in the ZFIM model, we will consider 
the gauge symmetry (or conservation of p(q)), but not 
require a constraint to be satisfied. 

Explicitly we can write x^r (or x' rr in the ZFIM 
model) in terms of the ladder sum, which is the solution 
to an integral equation (we define here various quantities 
to be used afterwards) 



J 



-y 



d 2 k 1 

r A(k, q,iu u )Q(k + -q, w n 



(2tt)= 

d 2 k . ,/(£k+4q)-/(£k-*q) 

2 A(k, q, UJ V ) — 



w„)£(k - -q,w„) 



(2n) 2 



Ck+iq £k-iq 



(4.26) 



Here A(k, q, iu) v ) is a one-particle-irreducible vertex function, 

1 . f d 2 k\ 1 1 

A(k,q,iw I/ ) = 1 - J ^^^(ki + -q,w n + w^)^(ki - -q,w„)r(ki,k,q,iw„) 



d 2 kx /(4 1 + i q ) - /(£ kl -i q ) wi . . > 
1_ ' 7^27 t r(ki,k,q,iov), 



(4.27) 



which we have written in terms of the particle-hole scattering series (the ladders with external Green's function lines 
removed) , 

T(k, k', q, iu v ) = V(k' - k) - i J |^ r ( k ' kl ' q ' *^)^( k i + 5* + Wv)a(ki - iq, u n )VQLi - k') 



V(k' - k) - 



d 2 k . . /(Ck+iJ-Z^k-lq) -^^ 

-— k.ki.q,!^ : — ^(ki - k . 

( 27r ) Ck+iq-?k-iq-^ 



(4.28) 



(Note that, in this approximation, the scattering func- 
tion depends only on the difference oj v of the Matsubara 
frequencies in the external fermion lines, and this is why 
we are able to perform the frequency sums explicitly.) 

Before analyzing these equations in detail, we pause 
to point out that for q, u> = iu u small and real, they 
have the form standard in Fermi liquid theory (see Pines 
and Nozieres (PN) f30|, and especially Nozieres fil) for 



the full, formal treatment), with the approximation that 
the /kk' function on the Fermi surface is taken to be 
t he lowest-order approximation as already given in Eq. 
(4.5), for spinless fermions, and this is just the content of 
the generalized HF approximation (see PN, Ch. 5). The 
Landau parameters Ft are then given by 



F e =Af(0)f e , 



(4.29) 



17 



ft = 



2tt 



/kk' cos £6<kk' 



(4.30) 



k F cos 6»kk' 



for 



for I > 0, where, as before, k • k' 
|k| = |k'| = kp. In particular, we notice that, since 
the density of states at the Fermi energy Af(0) = to* /2n, 
and since the bare kinetic energy is zero, comparison with 



Eq. (4.8) yields 



Ft 



This is a particular case of the relation 



m*/mo = 1 + Fi 



(4.31) 



(4.32) 



in ordinary two-dimensional Galilean-invariant Fermi liq- 
uids with bare mass m$. We can view the ZFIM model 
as such a system but with to — 00 ! from which F\ = — 1 
follows. This is the value that would usually be inter- 
preted as the borderline of stability of the system; how- 
ever, usually this view is taken because the bare mass 
is finite and the effective mass vanishes, and the latter 
causes instability. Here the effective mass is finite, so the 
system is not unstable, and moreover is held right at this 
point by this symmetry. We take it as implying that the 
ladder series must be analyzed with even greater atten- 
tion than usual to the limit lo — > 0, q — > 0, particularly 
for the 1=1 angular mode. We also point out a contrast 
with HLR, where this formula was invoked, but with the 



bare (or band) mass to in place of too, and was connected 
with Kohn's theorem and the f-sum rule. There the in- 
teresting limit was to — > (to send the cyclotron mode 
to infinite frequency), rather than oo. The present dis- 
cussion is clearly distinct, though it must be related at 
some deeper level. 

In Fermi liquid theory, relations like that above are de- 
rived through Ward identities connected with symmetries 
of the problem, and the symmetries are global, so the re- 
lations are most useful only at small q or w. Next we 
will derive a Ward-identity relationship between A and 
the self energy S within the HF approximation, in a way 
more directly connected with the symmetry generated by 
the p R 's, and valid for all w ^ and q. 

First we express the HF approximation as a pair of 
self-consistent equations: 



G(k,uj n ) = [iu n - (S(k) - A - /i)] _1 
1^ / d 2 k Y 



E(k) 



(2tt) s 



(4.33) 
V-(k-ki)0(ki,w n ) 



H^k-kO/^J 



(4.34) 



where £k = E(k) — /i — A as before (the direct term has 
been dropped as it plays no role in the following, for the 
one-interaction irreducible functions; it is absent anyway 
for the long-range interaction case). Then 



E(k+iq)-E(k-iq) 



x ( E(k a + -q) 



(V(k + iq - k x ) - V(k - iq - k!) ) c7i k, . ,•„ ! 

cPki ~ 1 
— — -V(k - ki)<?(ki + -q, u n + u v ) 

(Z7T) Z Z 



S(ki - -q) 



Q(ki - -q,uj n ), 



(4.35) 



after shifting dummy variables. But from Eqs. 
(4.27 4.28| ), — io; w A(k, q, iu> v ) obeys the same inhomoge- 
neous integral equation, the solution of which should be 
unique, so we conclude that 



iw„A(k, q, iuj„) 



,„-E(k+-q) 

Ju - Ck+iq + 4 



S(k--q) 



(4.36) 



which is the desired Ward identity (compare Ref. p8|). 
The left-hand side is the vertex function for dp R (q)/dT, 
which should vanish since p fi (q) commutes with the 
Hamiltonian. It implies that, if A is viewed as the scatter- 
ing amplitude for a fermion scattering off a potential cou- 
pling to p R , or for creating or destroying a particle-hole 
pair, then the amplitude vanishes if both particles are 
on shell, that is if their frequencies icu n satisfy icu n = £k- 
This suggests (following a similar argument in p2j , that 



was inspired by Q) that in the on-shcll states (energy 
eigenstates), if they satisfy the constraints G(q) = 0, 
then the latter property is actually preserved in the time 
evolution, in spite of its apparent violation in the HF 
states. This of course is because the calculation we have 
done is not the naive one of looking at the states as non- 
interacting particles, rather we used the conserving ap- 
proximation. It appears that the fermion excitations can 
after all be viewed as real physical excitations, satisfying 
the constraint conditions on physical states, even though 
the operators c* are not gauge-invariant and so would 
connect invariant to noninvariant states. These physical 
fermion excitations, which are dressed by the fluctuations 
around the HF states, are the physical composite or (as 
we s hall see) neutral fermions discussed in Jl2| and in 
Sec. ETBl 



Now we return to our original goal of calculating xWr 
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in the ladder approximation. Using the Ward identity 
and Eq. (4.26|), and assuming u ^ 0, we find 



LlUy 



d 2 k 

(27T) 



^2 (/(Ck+Iq) - /(Ck-iq)J 

(4.37) 



Another response function containing p fi that should 
vanish is XiSXl) ii*>„). In this case, the appearance of 
p L in place of one p R implies that the phase factors do 
not all cancel, and on using the Ward identity for the p R 
vertex we obtain 



d 2 k 



(2tt 



> l/(fk+iq) 



D ikAq 



0, 



(4.38) 



since shifting k by =F^q has no effect on the phase factor. 

As promised we have shown that the conserving ap- 
proximation guarantees that there are no fluctuations 
in p R (c[), at least for nonzero frequency. For zero fre- 
quency, the Lagrange multiplier fields A(q) (or the sub- 
set of diagonal elements, according to the final gauge- 
fixed form) enter to give the same result, but we will 
not show this explicitly. Similar issues were addressed 
extensively in the literature on slave bosons and heavy 
fermions in the 1980's (see for example ]6^ , |62] , |65[ -|68[ ) , 
and later in connection with theories of high T c supercon- 
ductors and quantum magnets. These problems also in- 



volve constraints, but these are usually abelian and gen- 
erate only U(l). It is still frequently stated incorrectly in 
the literature that in the functional-integral saddle-point 
approach to such problems, "the constraints are satisfied 
only on the average". In fact, as was well-known to sev- 
eral workers (such as the cited authors) in the field in the 
1980's, the correct RPA or 1/N (i.e., conserving) treat- 
ment of fluctuations yields just the same sort of results 
we have just derived, namely the vanishing of the vertex 
function for, and of all correlation functions containing, 
the constraint operators (like our G(q)), to all orders in 
the fluctuations. Thus the average of, and all fluctua- 
tions in, the constraints vanish, which means that the 
constraints are satisfied in every order of approximation, 
when this is set up correctly. (The extension to all orders 
for the present problem will be discussed later.) 

It remains to examine Xll- This will be undertaken 
in the next two subsections. 



D. Asymptotics of the ladder series 

In this subsection we continue the analysis of the con- 
serving approximation of the last subsection. We exam- 
ine the behavior of the ladder series at small q and w„, 
first to elucidate the mechanism behind the vanishing of 
X^r, an d then, in the following subsection, the results are 
applied to the calculation of the physical density-density 
response function xYl- 

The equation for T can be rewritten 



d 2 k x 



(27r) 2 <5(k'-ki) + f(k'-ki) 



/(£ki + iq) - /(fkx-iq) 



£k!-iq 



r(k,kiq s »a;„) =V(k-k') 



(4.39) 



which shows that it is a Fredholm integral equation, 
where the integral kernel appears in the curly brackets 
on the left-hand side, and contains q and lo v as param- 
eters. It implies that r is V times the inverse integral 
operator. The inverse could be calculated by finding the 
eigenvalues and eigenfunctions of the integral operator 
on the left. 

At ioj v — 0, (which could be viewed as the limit 
iu v — > 0), one zero eigenvector can be found for all q 
by use of the Ward identity proved in the p revious sub- 
section; it is £k+± q — Ck-J-q ( see Eq. ( 4.35 )). Thus for 
small iuj v , we expect to have, for all q, an eigenvector 
approximately £k+±q — £k-iqj with eigenvalue tending 
to zero with iuj v . If q — > also, we get 



4+iq-Ck-i q -q- v k 



(4.40) 



where Vk = Vk£k- At small q, the nontrivial part of the 
integral kernel becomes 



nk'-kx) 



de 



(4.41) 



which for zero temperature T is concentrated at k = kp 
(indeed, for all q, the difference of Fermi functions is non- 
zero only in a shell of width of order q around fcj? ). But 
this limit of the kernel is independent of q, so in addition 
to the eigenfunction just found which is proportional to 
cos 0k on the Fermi surface, there is another proportional 
to sin 6k. Note that these eigenfunctions, in the spirit of 
a Fermi-liquid analysis in terms of Sn^ or a deformation 
of the Fermi surface, are just rigid displacements of the 
Fermi sea, respectively parallel and perpendicular to q. 
The second eigenfunction is not a zero mode for q =/= 0, 
so is expected to acquire an eigenvalue that is nonzero as 
iui v — > 0, but vanishes as q — > 0. 

For general values of the ratio ito^/q the integral equa- 
tion and the eigenvalue problem are not easy to analyse, 
even for iui v , q small, where the eigenvalue equation takes 
the form 
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A(k,q, iu) v ) 



q- v kl 



of 



q • v kl iuj u de 



A(ki, q, = A(q, iu„)A(k, q, iw,,). 



(4.42) 



This form of equation is standard in Fermi liquid the- 
ory, with V(k — ki) replaced by — /kki- At T = 0, 
df /de — — <5(£k) and the equation can in principle be 
solved for k on the Fermi surface, and these values of 
the eigenfunction determine it elsewhere. Accordingly 
we might expand both A and V in terms of Fourier 
modes cos£6>k, sin^# k , £ = 0, 1, . . . , for |k| = kp. For 
iui v /\<\\vp ^ 0, the Fourier modes are mixed by the inte- 
gral kernel, so that all components of 

oo 

- V(k - k') = fo + 2 ft cos^kk* 

oo 



= fo + 2 M cos ^ cos 



1=1 



+ sin^ k sin Ifa) 



(4.43) 



are involved. We have seen that the 1=1 mode and f\ 
are crucial to the analysis and must be kept. The other 
Landau parameters F e take no special values and merely 
produce finite renormalizations of the response functions 
(some identities are implied by the existence of the zero 
mode for all q, but these bring in derivatives of Vk and 
thus parameters that lie outside of Fermi liquid theory) . 
We propose just to drop these effects so as to obtain the 
simplest possible approximation that is still conserving. 
This can be done by replacing fi for I ^= 1 by zero, or 
more accurately by assuming that the only eigenfunctions 
A that are needed are just q • v k /g, q A v^/q (which are 
the correct continuations off |k| = kp). We will actu- 
ally use this even to higher order in q, as we will see is 
necessary. 

With this further approximation, the eigenvalues cor- 
responding to the two eigenfunctions can be evaluated. 
The final result for T is 



r(k, k', q, iu v ) 



q • v k q • v k ' q A v k q A v k ' 



uixo (q, iuv) q 2 xo (q> *w„) 



, (4.44) 



where 



Xo(q,M^) 



d 2 k /(4+iq)-/(4-Iq) 



( 27r ) 2 fk+iq-fk-iq 



(4.45) 



is the "density-density" response function of a Fermi gas 
with dispersion £k, and 



Xo (q, Wu, 



1 Knn\ 2 f d2k f <l Av A ' 



d 2 k 

(27T) 2 

/(C k+ i q )-/(^-iq) 



(4.46) 



is the transverse "current-current" response function of 
the same Fermi gas, including the q-, ia^-independent 
contact ( "diamagnetic" ) term, xo arose in a similar way 
from the longitudinal current-current response, on using 
the continuity equation. Note that what we are calling 
the "density" and "current" , though natural in appear- 
ance, are not to be identified with the physical density 
and current. 

The above expressions for \o and Xo" are valid for any 
q and iuo v . On the real frequency axis, at ui/qvF and q 
small, they become 



Xo (q,^ + i0+) = Af(Q) + iAf{0)uj/(qv F ) 



Xo (q. 



i0 i 



q 2 X*d + iuk F /(2nq). 



(4.47) 
(4.48) 



Here Xd is the diamagnetic susceptibility of the Fermi 
gas with dispersion £k- It is a non- Fermi- liquid prop- 
erty that involves derivatives of Vk at kp; if £k were 
= (k 2 - fc|)/2m*, then x*a would be = -l/(127rm*). 
These imply that the eigenvalues of the longitudinal and 
transverse eigenmodes of the integral kernel above vanish 
in the ways predicted in this limit. This involved the can- 
cellation of the diamagnetic term in the current-current 
response in both cases; this cancellation is well-known in 
normal fluids (i.e., non-superfluids). 

We can now show that even this further a pprox imation 
is conserving in the sense discussed in Sec. 1VC. Using 
the above form of T, we can calculate 



X 1 rr 



Xo - Xo(xo) 1 Xo = 0, 



(4.49) 



where the second term is the contribution of T, for all q 
and iuj v ^ 0. In this calculation, the transverse mode in 
r did not contribute. A similar calculation shows that 
Xrl = 0- An exact treatment of the ladder series in the 
regime uj/qvp <C 1 and q <C kp yields the same form with 
all xo's replaced by Xo/(l + Fo), and the cancellation still 
occurs, in agreement with the previous subsection. 
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E. Physical response functions 



1. Density-density response function 



In this subsection we calculate Xll> tne physical 
density-density response function, and its limits, the 
compressibility and longitudinal conductivity. We also 
consider the scattering of the fermions by an external 
potential, and the expression for the current density. 



As already remarked, the fact that the p L vertex con- 
tains the opposite phase factor from that in p R means 
that not all the phase factors cancel in Xllj instead those 
at the two vertices at the ends of the ladder are doubled. 
We have 



irr f dkdk /gkjjqj /gk-jq) , , . /gg+jq) /gg-frq) ik A q _ ik ' Aq 

Xll = Xo + / ,„ x 4 7 ; : — 1 (k, k , q, iw„J - ; — e 

However, by comparison with Xbm = X'rl — Xl"r = 0> this simplifies to 



Xll 



^ ( e 'ikAq _ jWg-ikAq _ j\/(£ k +iq) /(£k-Aq) 



(4.50) 



(2tt) 2 
xr(k,k',q,«w I/ ) 
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jkAq 1 /(^ + lq)-/(ek-lq) 

(2tt) 4 1 J ' 



Ck+iq £k-iq 



lU>v 



/(4' + ^q) A^-|q) , _ ik / Aq _ 1 - 
4' + |q-Ck'-iq-^ 



(4.51) 



For small q, we now expand the phase factor. The first 
term is then the form found by [^-^6| . It is the same as 
putting p L — p R in place of p L , which goes as ~ ik A q 
at small nonzero q. The second term is the ladder series 
with the insertion (kAq)(k' Aq) at the two vertices. This 
exhibits the effectively dipolar nature of the coupling of 
an external scalar potential to the physical density: the 
fermions carry a dipole moment Ak, as found in Refs. 
2|,f24h26 and discussed in Sec. [ll| In V, only the trans- 



verse mode now contributes, and we obtain 
Xll = q 2 m*(p + m*Xo(q, iuv)) 

Xo"(q, iui v ) 

= -q 2 p{p + m*Xo(ci,^))/xo(<l,^). (4.52) 

Note that in the numerator, the p's occur because of 
the absence of a "diamagnetic" term to cancel it, and in 
writing the remainder of the numerator as Xo we have 
neglected the difference between k/m* and v k , which af- 
fects the coefficient of the term in Xo quadratic in q. 
This term can be neglected anyway in the following. In 
the small lo j (qvp), q region we then have 



~X*d - iuk F /(2irq 3 )' 



(4.53) 



This is similar in form to the result obtained by HLR, or 
the renormalized version of it according to the scenario 
discussed in Sec. II A, if we note that p = 1/(2^0) in 
general (and </. 



I A, if we note that p 
1 here), except that the 1 in the de- 
nominator in Eq. (2.5) has been dropped. That 1 came 
from the Chern-Simons term, which couples longitudinal 
and transverse fluctuations; by contrast, in the conserv- 
ing approximation in the present approach, the ladder 



propagator T does not couple these modes. Note that 
the first term in the first line of Eq. ( 4.52| ) is essentially 
the result of Refs. &p6|, 



irr 

Xll 



q 2 m* (/> + m*Xo"(q,iw„)) , 



(4.54) 



which behaves differently at low lo and q, as we will see. 

We now take various limits of this expression. As 
lo — » 0, we obtain 



^=iim x I £i(q,o) = 

dp |q|—0 



5 2 /x3, 



(4.55) 



which is finite and positive, so the system is compress- 
ible as in HLR, though agai n the expression differ s from 
that in the scenario of Sec. [I A , as given in Eq. ( 2.13 ). 
Though we used the approximate form for L, our result 
is exact within the ladder (conserving) approximation. 

To obtain the low-frequency longitudinal conductivity, 
relevant to the surface acoustic wave experiments, we de- 
fine a relevant limit : 



Czx(q) = lim 

ui/q— >0 



lim 

uj I q fixed 



-X^(q^ + *0 + ) (4.56) 



for q parallel to x (the conductivity should always be 
viewed as the response to the total electric field, so it is 
related to the irreducible response). Here "lim" means 
that we keep the leading nonzero term. This limit cor- 
responds to considering a long- wavelength sound wave, 
so |q| is small <C &f and lo = |q|f s , and then taking the 
sound velocity v s to zero, (i.e. v s <C vf)- Then we obtain 



,2irq 



2nk ? 



(4.57) 



in exact agreement with HLR for <j) — 1. There a differ- 
ent procedure was used to define a xx {q), as given by HLR 
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eq. (B4.a). That and the present definition give the same 
result both in the RPA of HLR and in the present ap- 
proximation. This result was expected to be very robust 
on Fermi liquid grounds, within the scenario discussed in 



Sec. [I A , since it corresponds to the transverse conduc- 



tivity of an ordinary Fermi liquid, which is unrenormal- 
ized in Fermi liquid theory. Remarkably, it is the same 
here, in spite of other differences in the structure of the 
expressions. This result is not obtained from the expres- 
sion (4.54) p4| . It is also remarkable how the factor p, 
which came from a standard gauge-invariance result for 
the usual Fermi liquid, here plays one of the roles played 
in the CS theory by a xy (= p in our units). This effect, 
that the "current" response at oj/q — > of a Fermi gas to 
a scalar potential coupled to the dipolar expression for 
the density gives the Hall conductivity, was pointed out 
by Stormer p9| . 

Fi nally the spectral density for x'£i(q,w) implied by 
Eq. (4.53) is at low frequency 



xTl (q,«) 



ujk F p 2 /(2irq 3 ) 



3\2 



(4.58) 



(but vanishes for \uj\/{c[Vf) > 1) and has a peak, an over- 
damped mode at ui ~ |q| 3 , similar to the result of HLR. 
As many physicists have noticed, this implies for the var- 
ious moments, as q — ► 0, 



oj n ~ q n+3 , n > 1; 

~ q 3 lnl/g, n = 0; 
~ const n = — 1. 



(4.59) 



For n < — 1, the moments diverge as usual. 

The n = 1 moment can be obtained exactly, because 
of the Kramers-Kronig relation, 



— OO 

1 

2 



du< xT L (qV) 
7T u/ - (oj + i0+) 



dJ 



u'^Il (qV), (4-60) 



as oj — > oo. The high-frequency behavior of &t 
small g can be obtained by returning to the integral 
equation for T, Eq. ( |f.39] ). To leading order in qvp/ui, 
r(k, k', 0, oj) = V(k — k'), and we obtain from Eqs. 
fl4.43|j4.5]f) 



fk 2 F p(l 



Fi 



4oj 2 m* 



(4.61) 



(The same result except that F-i is replaced by zero is 
obtained using our earlier approximation for T.) This 
can be compared with the result in a usual Fermi liq- 
uid, which is —q 2 p(l + F±)/ '(oj 2 m*) — —q 2 p/(u 2 m) on 
using 1 + F\ = m* I'm. We return in Sec. [v| below to the 
question of the general validity of our result, beyond the 
ladder approximation. 

The moments of the spectral density of the full re- 
sponse function xll can now be obtained also. For the 
n = — 1 moment, one finds ~ V^q) -1 for a long range in- 
teraction, as usual in a compressible system. The n = 
moment behaves as g 3 lnl/g again, and gives the LLL 
"static" (equal time) structure factor s(q). It does not 
go as q 4 as GMP suggested it should in any liquid state. 
This is because compressible liquids have both low-energy 
modes and long-range correlations that produce nonan- 
alytic behavior of s(q). GMP concluded that fluids in 
the LLL should be incompressible, but this argument is 
invalid (this point has also been made by Haldane |f27f). 
The n = 1 moment goes as q 4 , as argued by GMP, and 
using the high frequency behavior of XLi,(q, w), an d be- 
cause V(q) is less singular than q~ 4 , 



du' i „, f°° duj' ■ it 

—uj xll (q, w ) = / — w x LL (q, w ) 



oo 
4^,2 



q 4 k 2 F p(l + F 2 ) 

8m* 



(4.62) 



to leading order in q. GMP found a formula for this mo- 
ment in terms of V(q) and s(q), so we obtain a relation 
among the quantities to*, F%, and s(q). The result for 
the n = 1 moment of can also be viewed as a sum 
rule for the leading part at small q of the longitudinal 
conductivity Re a xx (q, oj) = ojXll" w )/<? 2 - 



2. Fermion scattering vertex 

We now consider the scattering of the fermions by an 
external potential V^ x t(r,t). The scattering of a fermion 
from wavevector k+iqtok— iqis given in the same 
ladder diagram approximation by the vertex function, 
similar to A earlier except for a phase factor, 



A L (k,q,^) = e lkAq 



d 2 k x 



5 ikiAq 



/(Ckj + iq) - /(^ki-iq) 



r(ki,k, q, iw,,) 



(4.63) 



after removing the same phase on the external lines as 
for A (only the irreducible part is shown). If the phase 
factors are replaced by 1, we obtain A, so we will first 
reconsider this briefly. 



Earlier we showed that, in the small q limit, 

A = 1 - q • v k /(iw„). (4.64) 
In terms of the asymptotics of T, the second term is the 
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correction produced by the longitudinal mode. While 
the first term is the bare scalar coupling to the exter- 
nal potential, the second term couples to the fermions 
through their velocity, that is to the "current" (in the 
same sense as before), and so can be viewed as describ- 
ing a longitudinal vector potential. Because of the factor 
q/iuj v , the vector potential cancels the direct effect of the 
scalar potential, if we consider the electric field they pro- 
duce. The system responds by producing a longitudinal 
response purely in the form of a vector potential, because 
we chose the gauge such that the scalar potential in the 
functional integral vanishes at nonzero frequencies. Thus 
for gauge-invariant response functions, such as XRR that 
we considered earlier, these terms produce complete can- 
cellation, as we saw earlier in the example. This should 
also be true in other calculations, such as for the effect of 
an external "impurity" potential on the conductivity, if 
it coupled to p R instead of to p L as it would in fact (such 
an "impurity" potential would be static, but as usual 
the same effects would be found there as for all nonzero 
frequencies, thanks to the zero-frequency Lagrange mul- 
tiplier or scalar potential field). 

Since the vertex functions A and A L differ only by 
the phase factor, we conclude that the phase factors like 

e «kAq can kg r gpl acec J by g«kAq _ ^ wnen using A L . To 

first order in q, this gives the dipolar coupling k A q with 
dipole moment Ak. The first term in A L is thus the direct 
coupling of V^ x t to the dipole moment of the fermions. It 
should be contrasted with the direct, minimal coupling to 
the fermions with charge 1 in the scenario for the low en- 
ergy behavior in the approach of HLR , described in Sec. 



II A. In the second term in A L , where the ladder series T 



contributes, the dipolar coupling brings in the transverse 
mode in the ladder series, as in the calculation of Xll- 
This coupling gives essentially 



q A v k ( i o + m*Xo(q, iWi/))/Xo (q, iu) v 



(4.65) 



at small q, iui v , which is a coupling to the transverse 
current, and is similar to that found in HLR and also in 
|jlf in connection with the effects of an impurity po- 
tential, that is the iu) v = limit. As there, the ex- 
ternal potential couples to the density, which induces 
a transverse vector potential, which, because it is sin- 
gular at q = 0, scatters the fermions much more effec- 
tively than the direct minimal coupling to the potential, 
let alone the dipolar coupling. The scattering produced 
can be simplified by comparison with the physical den- 
sity p L induced by the same external potential, which is 
(P L ) - P — XLiKxt(q, iu>u)e~i^ 2 . This shows that if 
the induced transverse vector potential is denoted a + A, 
then we have 



VAa 



<P L )/P 



-2-kc 



(4.66) 



which is exactly the equation in the CS theory! This 
shows that the fermions experience a vector potential that 



obeys eq. ( 4-6t ), where p L is the physical charge density, 
even though there is no CS term in the effective gauge 
field coupling and the fermions behave as dipolcs. This 
agrees with the use in Ref s. Jl0| , ^2| of the Berry phase 
argument of Arovas et al. ]23|] to obtain the vector po- 
tential seen by the fermions, which in no way assumed 
that there are flux tubes attached to the particles, unlike 
the CS approach. Note that, since we also have 



P 



P - V A g, 



(4.67) 



this is consistent with a + A = g/p for the longitudi- 
nal part. There should also be an equation —a — Vao = 
2-K(f> A j L , where j L is the physical current density. The 
problem of the form of j L in the present approach will be 



considered in Subsec. IV E 4 



3. Effect of impurities 

Here we consider the effect of impurity scattering on 
the density-density response and the longitudinal con- 
ductivity. The HF and ladder approximations can be 
reconsidered with impurities present. We neglect here 
the mechanism of the preceding subsection, and take 
only direct scattering by the impurities, analogously to 
the bare HF considered so far. The average self energy 
should contain an impurity line (the self-consistent Born 
approximation), and the ladders contain both impurity 
lines and interactions as the rungs of the ladder. The ef- 
fective mass and the diamagnetic susceptibility will gen- 
erally be renormalized by the impurity effects, but we 
will not distinguish them from their counterparts in the 
pure system. Calculations are straightforward, and the 
results can be written down using well-known formulas. 
The scattering rate 1 /r is given by the usual expression, 
but contains m* from the density of states (this could be 
replaced by the rate from the mechanism of the preceding 
subsection, but this makes little difference). At q = 0, 
we have 



iwp{p + m*Xo) 



Xo 



(4.68) 



and, in the Drude approximation, recalling that the 
current-current response is isotropic at q = 0, 



X^(0,^ + z0+) 



iLupr 



m*(l — iuor) 



Then 



c (0,lj) = pm*/r = cr 



(4.69) 



(4.70) 



independent of to. This can be viewed as the usual form 
of resistivity of the fermions, p xx = (pr/rn*)^ 1 , divided 
by p xy , so is consistent for small p xx with the result of 
the CS theory, of adding the fermion and CS resistivities 
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(see Eq. (2.9)). The frequency- independence is also con- 
sistent with this, if in the CS approach one uses m* in 
place of m, and includes FL corrections as in the scenario 



described in Sec. [I A. The effect of the latter corrections 



is to replace 1 — iujr by 1 — iujTm/m* (see PN, p. 191). 
As m/m* — ► 0, with m* , r fixed, the result above is 
obtained. 

For finite wavevector, we will consider only the small 
uj and q region. With impurities present, Xo is analytic 
in q 2 and uj, 

Xo{<l,u + iO + ) = q 2 X* d + iupT/m*. (4.71) 
We then obtain the longitudinal conductivity 

iujao 



a xx (q,u + = 



(4.72) 



iui — Dq 2 ' 

which exhibits a diffusion pole, with diffusion constant 

D = -m* X * d /(pr), (4.73) 
and Co obeys the Einstein relation <tq = Ddn/dp. 

4- Physical current density 

We turn here to the calculation of the expression for 
the physical current density within linear response. The 
most obvious way to obtain the current is by projecting 
the usual expression to the LLL, as was considered by 
GMP. This yields 



AVp L /(2m), 



(4.74) 



which involves the bare mass, and describes the current 
due to the cyclotron motion of the particles. Since it 
clearly obeys V ■ j c = 0, and gives zero when integrated 
across a section with a boundary condition of zero den- 
sity, it does not contribute to transport. This current, 
when coupled linearly to a change in the vector poten- 
tial, A-j c , describes a magnetic moment on each particle, 
which should be recovered in the U(l) CS approach, as 
argued by the authors of Ref . , and obtained by SM 
& 

We are concerned with transport and with response 
functions, and this part of the current contains explicit 
derivatives, so is of less interest at long wavelengths. We 
therefore turn to the current due to drift motion of the 
guiding centers of the cyclotron orbits of the particles, 
due to both the external one-body potential Aq, and the 
interparticle two-body interaction. We will not consider 
fully the response to a change in the physical vector po- 
tential A. The existence of both parts of the drift cur- 
rent was recognized by GMP and in Ref. [JlO|; for fur- 
ther discussion, see Refs. ]7l| , [72| . In principle, they can 
be obtained by carrying the calculation of the projected 
current to higher order in l/ui c (the cyclotron current j c 



being the leading term, of order w c ), by considering vir- 
tual excitation of the particles to higher Landau levels. 
This is carried out in Ref. |Q ; it yields two types of terms 
of order w° in the matrix elements of the current within 
the LLL, for an external potential V CX f The first of these, 
called j^, can be written as a series of derivatives of the 
LLL-projected potential V cxt and of the density p L ; the 
series can be further divided into a series of exponential 
form that agrees with the "Noether current" of Martinez 
and Stone, and another series, beginning with a third- 
order derivative, that is of the form of an integral of an 
exponential. The second type of term |7j| consists of the 
modification of the cyclotron current by the effective LLL 
Hamiltonian to order u~ , so is more complicated. The 
general expression for the current is thus by no means 
simple. However, to find the net current for transport 
purposes, we require only the small q limit, and for this 
the result is just 



j L = -P L A VA 



(4.75) 



for a slowly varying potential Aq — V cx t, which exhibits 
the Hall conductivity a xy — p in our system. 

For the small q drift current due to the interaction, we 
have in Fourier space 



j L (q) = 



d 2 q' 



:*AqT(q'):/(q + q')/(-q') 



(4.76) 



Diagrammatically, one can see that to calculate the lin- 
ear response current to a scalar perturbation within the 
conserving approximation, it will be sufficient to take the 
operator itself in the HF approximation. Since = 
in the unperturbed ground state, the leading term is ob- 
tained by replacing a pair of operators c\ c by their ex- 
pectation value in the ground state, 



(4 c k2 ) = (2 7 r) 2 <5(k 1 - k 2 )0(k F - ki) 



(4.77) 



in all possible ways, that is two "direct" and two "ex- 
change" terms. Of the direct terms, one vanishes and 
the other is seen to give the Hall current produced by 
the field due to the interaction with the average density 
of particles at wavevector q, 



j i (q)dircct = ~i AqpV (q)p L (q). 



(4.78) 



In calculating the irreducible response to the total field, 
this term is clearly included automatically. Therefore we 
can turn to the exchange terms which alone give the ir- 
reducible response. Since q is small, we use 

0(fc F -|k+iq|)-0(fc F -|k-iq|) 

= -qcos6 k 5(k- k F ) (4.79) 
for q in the x direction, and after some algebra we obtain 
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j L (q) iri = / A (k - k>(k - k')ct _ , , 

x (—q cos 9k' 5(k' — kp)) 
d 2 k 



k+iq 



[iAkq-k(l + F 2 )/m* 



+ i A qk 2 F (F - F a )/(2m*)] C t_i q c k+ i q , 

(4.80) 

where the Landau parameters Fg were defined earlier, in 
Eq. (4.43). We assumed that only values of k near kp will 
be used, which is true for linear response (thus k F = k 2 ). 
Interpreting cj t _i q c k+ ^ q as 8rik(q) in FL theory, 

where 5n k (r) is the departure of the distribution of occu- 
pied k values at r from the ground state, and is assumed 
to be nonzero only for k near kp, this can be identified 
as 



i£(q)irr = -ie^qxILuxiq) 



(4.81) 



where 

n„ 



(kph,- -k 2 8^){\ + F 2 )/m* 



d 2 k 
(2tt) 2 



<5n k (q) 



(4.82) 



is the stress or momentum flux tensor of the FL; it is 
equivalent to that in Ref. j74|, modified to two dimen- 
sions. Since we have identified p L (r) — p — V • P and 
P(r) = Ag(r), we expect a term in the current j£ r = P(r) 
p5[. But by momentum conservation, 



d 9^ , « rr 



0, 



(4.83) 



and so we find Eq. (4.81). Since we also wish to identify 
a + A = g/p, we find 



Jirr 



p A a, 



(4.84) 



which is essentially the other CS-like equation. 

We should also add to the Hamiltonian the potential 
terms 



d 2 q 



ao(q)p H {~q)+A p L (-q) 



(4.85) 



where do is the scalar potential introduced earlier, which 
implements the constraint p R — p, and for which we 
chose the gauge do = 0, and A n (q) = e~i q A (q) = 
Vcxt(q) is the externally applied potential. Then the 
right hand side of the momentum conservation equation 
becomes 



(pVa +p L (r)VAo) 



(4.86) 



at long wavelengths. Here the coefficient p arises from 
p R on using the constraint. There is also a similar Hall 



contribution to — p L AVAo to the current density j L . Ex- 
pressing the total physical current j L in terms of g = a//5, 
we obtain 



p A (a + Vo ), 



(4.87) 



which is manifestly gauge-invariant and of the CS form. 

If we consider the current in the right coordinates, j R , 
in a similar way, we find that in the absence of clq it 
vanishes identically, because p R commutes with H . This 
result of vanishing current was already invoked in Sec. 
II D. It can be interpreted by breaking the current into 
the pieces g/m* and (a + A)p/m* shown there. The 
first term represents the velocity of the fermions, while 
the second represents the usual backflow correction in a 
FL, which in the present case of F\ — —1 exactly can- 
cels the first part. The same effect occurs in the ZFIM 
model: the total current carried by each fermion is k/m 
by Galilean invariance, and m = oo, so it vanishes. (In 
the presence of do, we find j R = p A Vao, a Hall c urrent. 
This does not affect our argument in Sec. IIP, which 



uses only the irreducible part of the current, from inter- 
actions.) A similar calculation can be given for j L . The 
velocity term and the leading part of the backflow are the 
same as for j , and so cancel. The subleading terms then 
give the result as calculated above. This cancellation of 
the leading terms is (perhaps not surprisingly) similar to 
what occurred in the formula for the density p L on using 
the constraint on p R . 

The irreducible longitudinal current density-density re- 
sponse function xflpL should be ui/q times xYl- This 
can be verified in terms of the ladder series expressions 
for both, if one consistently either keeps or drops the 
Landau parameters Fg for I ^ 1 in both the ladder series 
and the q-j L vertex. In particular, in the small q/oj limit, 
the (1+ F2) J m* term in j L reproduces that in Xll- How- 
ever, if we consider the longitudinal current-current re- 
sponse, which should be uj 2 jq 2 times Xll,' we see ^ na ^ * ne 
two-point current correlation function starts at higher or- 
der in q/uj than the required term (two-point correlation 
functions always vanish as uj — > 00). A similar difficulty 
is familiar in the usual Fermi liquid and is resolved by 
the presence of a term in the current, —pA/m (the "dia- 
magnetic current"), that is linear in the applied vector 
potential perturbation, so that the response function (xg 
in the noninteracting case) consists of a constant plus the 
two-point function of the current without the A term. A 
similar effect should occur here. The term required in j L 
is of order q 2 . One might attempt to find such possible 
terms by making the stress tensor expression Eq. (4.82) 
gauge invariant by replacing all k's (including k 2 F = k 2 ) 
by k — a — A. This does not affect the other calculations 
done up to now because, in the absence of a perturbation 
in the external A, the net a + A does not contribute in 
linear response. But further work is required to check 
the form of this tensor, since the gauge invariance under 
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SU(iV) or Woo reduces to conventional U(l) gauge the- 
ory only at long wavelengths, while this expression for j L 
is higher order in derivatives. In any case, such minimal 
coupling terms do not produce the necessary factors of 
q, and so should be absent. A way to find the part of 
the longitudinal current linear in a change in A, which 
should be correct at long wavelengths, is to add a term 



-5A ■ (where j i< -°- > is the exact expression (4.76) 

of zeroth order in the perturbation <5A) to the Hamil- 
tonian, then calculate the longitudinal current through 
first order terms in 8A by commuting p L with H . The 
resulting first order term can be seen to give the cor- 
rect high-frequency limit of the response, because it is 
given by a double commutator of H with p L (±q), which 
is what appears in the sum rule for the first moment of 
the spectral density of Xllj an d we have seen that it is 
related to (1 + F 2 )/m* also. Thus the correct term is 
obtained, and must be used in the longitudinal current- 
current response for all co/q to ensure agreement with the 
density-density response. 

We now consider the full conductivity tensor at q = 0. 
The longitudinal part has already been considered. The 
full conductivity tensor can be written in the Kubo form 



a^(0,uj + i0 + ) = pe^ + 



i{uj + i0+) 



xg^(o,a, + ;o + ) 

(4.88) 



where the first term is the Hall conductivity and x^lal is 

J pi J V 

the current-current two-point function for the irreducible 
part of the current. This form was proposed by Lee [fj5| . 
We may also consider the conductivity tensor when im- 
purities are present. Note that the q 2 term in j L does 
not contribute when q = 0, even when impurities are 
present. However, we expect an additional contribution 
to j L from the impurity potential, which we have not 
explicitly calculated. Because averaging (using Gaussian 
disorder) produces diagrams like those for interactions, 
except that no frequency is transferred along impurity 
lines, it should be similar to that derived above. It will 
represent the loss of conservation of momentum when 
disorder is present. Only the off-diagonal part of X^l^l, 
or the corresponding transverse response to a scalar per- 
turbation, has not so far been calculated. Because the 
ladder diagrams in the interaction and impurity lines do 
not violate parity (reflection symmetry), there can be no 
off-diagonal terms unless the impurity current vertices 
that we have not calculated contain pieces both parallel 
and perpendicular to q. If such terms are absent, then 
o~xy = P, unaffected by impurities in this approximation. 
As emphasized by Lee p5[, this differs from the result 



of the U(l) CS approach mentioned in Sec. [I A. It was 
argued in Ref. [j75| that in the U(l) CS fermion approach, 
applied to the v = 1/2 case, particle-hole symmetry im- 
plies that fj xy — 1/2 exactly, which is only satisfied by the 
scenario described in Sec. II A if a^, xy of the CS fermions 



is —1/2. Assuming our results also apply to v = 1/2, 
there is clearly no problem with particle-hole symmetry 
in our self-consistent Born approximation (SCBA). We 
should point out, however, that in this or the similar 
approximation for the U(l) CS approach, the results do 
agree at leading order in p^, xx / p xy , and the condition 
o~ip X y = —1/2 is only needed to guarantee a xy = 1/2 
to all orders in this expansion. Thus the contrast be- 
tween the naive SCBA result a^ xy — and the required 
o~ip xy — —1/2 is not such a dramatic singular correction 
as it might appear at first sight. At higher orders there 
will of course be other correction terms not included in 
the SCBA, which can drive the system into the criti- 
cal regime representing the transition between quantized 
Hall plateaus. 



V. EXTENSION TO ALL ORDERS IN THE 
INTERACTION, AND DISCUSSION 

In this Section, we consider the extension of the results 



of Sec. IV to all orders in the interaction, and describe 



the structure of the results we expect, in a scenario which 
replaces the previous U(l) CS scenario described in Sec. 
[I A. First we consider a more complicated conserving 



self-consistent approximation, with special attention to 
long-range interactions. Then we explain the FL theory 
structure for sufficiently long-range interactions. 

In the HF and generalized HF approximation of Sec. 
[V , the exchange diagrams contained the bare interaction 
V(q), and this led to a vanishing m* at fcp for Coulomb 
or longer-range interactions. An obvious improvement to 
make is to ins ert the ladder series into the Coulomb ver- 
tex, as in Sec. [V E 2| . The longitudinal part of the ladder 
series T renders the coupling to the fermions dipolar at 
long wavelengths, which removes the divergence in 1/m* 
for interactions less singular than l/q 3 . At the same time, 
we can insert the ladder series inside the interaction line 
itself, thus screening the interaction. We can also re- 
place the interaction line in the exchange diagram by T. 
Finally, we make this approximation self-consistent by 
making these replacements for all interaction lines, in- 
cluding those in T, thus iterating to self-consistency. This 
approximation, applied to response functions as well as 
the self energy, is once again conserving in the same sense 
as in Sec. IV C , and the conclusions there, which follow 



from Fi = —1, still apply. 

This approximation is clearly not as tractable as HF, 
but we can still make some general statements. The sys- 
tem should still be compressible for all interactions con- 
sidered (those less singular than 1/q 2 as q — > 0). The 
longitudinal mode in the ladder just produces the dipolar 
coupling effects already mentioned, which do not cause 
a breakdown of FL theory, though the effect of the ex- 
change self energy that contains T in place of V has not 
been calculated. The transverse mode in T produces 



2G 



singularities in the self energy for Coulomb or shorter- 
range interactions. The self-consistent summation pro- 
posed here is the same as regards the transverse mode as 
that studied in Refs. fll]JH] (and similar to that in Ref. 
p7[ ) . We have nothing to add here to the previous discus- 
sion of this case, except to emphasize that these singular 
effects should be treated after the other FL renormaliza- 
tions discussed in this paper, and that in relation to the 
U(l) CS approach, the effects incorporated in this paper 
are related to the longitudinal, not transverse, CS gauge 



field fluctuations (see Sec. II A). For interactions longer 



range than Coulomb, there is no breakdown of FL theory, 
since to* remains finite and the quasiparticle decay rate 
vanishes faster than the renormalized excitation energy 
£k as k — » fci?, though not as fast as k 2 . 

We can now discuss the general structure expected in 
the results to all orders in the interaction; some of this 
is implicit in the foregoing discussion. We consider only 
interactions longer-range than Coulomb, so there is no 
breakdown of FL theory. For Coulomb interaction, the 
results are probably still useful, since the only other effect 
is a logarithmic divergence in to*, which is very weak. 

For such interactions, we again separate in the response 
functions the "direct" or reducible diagrams, which rep- 
resent the long-range self-consistent field produced by the 
expectation value of the density. The remaining diagrams 
are analyzed in terms of the fermion-hole irreducible scat- 
tering vertex, which at q, u> — ► is nonsingular and 
defines the parameters fi and hence Fg = m*ft/2'K. 
A Ward identity, now valid to all orders, implies that 
F\ = — 1. In fact an identity for the p R vertex, like 



that in Sec. IV C, is valid to all orders and for all q and 
w„ ^ 0, and expresses the fact that [p R , H] = 0. In the 
general diagrams that contribute to these vertex func- 
tions, the phase factors in the interaction vertex do not 
all cancel, so the system is not equivalent to the ZFIM 
model. The results nonetheless have the same structure 
as in Sec. |i"v| , and at long wavelengths can be interpreted 
in terms of an infinitcly-strongly coupled gauge field, cou- 
pled to the FL. There are no parity violating effects in 
the long-wavelength dynamics of this system, because the 
Landau interaction /kk' is even under exchange of k and 
k' . The only parity- violating effects come in the coupling 
to external electromagnetic fields, where the Hall effect 
appears, and the physical density and current obey the 
CS-like equations. The self-consistent field produced by 
the long-range interaction (the reducible terms) also pro- 
duces Hall currents, but there is no parity violation be- 
cause interactions within the system couple to the density 
at both ends. The fluctuations in the longitudinal part 
of the gauge field can be reconsidered by changing to the 
gauge V • a = 0, in which it is the scalar potential ao that 
fluctuates (at all frequencies). This absorbs the Fq we 
had previously, and the condition p R = p is maintained 
through an effective Fq that is now infinite (the Landau 
parametrization is not gauge invariant). The longitudi- 



nal part of the ladder series at low ui/q gives an effective 
interaction between the fermions, which is of order the 
inverse density of states (this is similar to effects in the 
local Fermi liquid in the Kondo problem, see |M). Be- 
cause the leading "monopolar" part of the p~density 
fluctuations is suppressed by this, the leading nontriv- 
ial part is described by the subleading, dipolar part of 
the exact density expression p L (note that this sublead- 
ing coupling is not described by the minimally-coupled 
long- wavelength Hamiltonian Eq. ( 2.30| )). A notewor- 
thy feature of our approach is that this is not obtained 
separately from the transverse gauge field effects, nor in- 
serted at the beginning, but emerges later. The dipole 
moment Ak on each fermion is not renormalized, because 
the momentum is a conserved quantity. This really de- 
serves an explicit proof, but it will be omitted because of 
the similarity to results in standard FL theory (see, e.g., 
Nozieres f6l]| ); quite generally, conserved quantities are 
not renormalized. 

The compressibility is given by 



dn 
dp 



P_ 

X* d 



(5.1) 



where \d is the fully renormalized (irreducible) diamag- 
netic susceptibility, and is the only non-Fermi surface 
quantity to make an appearance in the response in the 
regime q, to small. The other quantities mentioned in 
Sec. IV are given by the same forms as there, when 
written in terms of p, kp, to*, Fg, and Xd- I n partic- 
ular, we mention the longitudinal conductivity in the 
regime q 3 < ui < qvF, relevant to surface acoustic waves. 
The result, which is identical to that of HLR, is ex- 
act in the same way, and for the same reason, as the 
low- frequency transverse conductivity of the usual FL. 
Also, the high frequency behavior, or n = 1 moment 
of the spectral density, of the irreducible density-density 
resp onse, is given by the same sum-rule-like form as in 
Sec. IV E 1 , as long as we consider only excitation of a 
single quasiparticle-quasihole pair (in the FL sense). If 
multiple quasiparticle-hole pairs do not contribute at this 
order in q, then this "sum rule" is exact. In the usual 
FL, multiple quasiparticle-hole pairs contribute to spec- 
tral densities at 0(q 4 ), by considerations of phase space, 
and the f-sum rule is for the q 2 part (and higher-order 
terms actually vanish in this particular case). Thus it 
is not certain in our case that our sum rule is exact. 
The same phase-space considerations apply, and if we 
assume that the squared matrix element of the density 
p L is of order q 2 (i.e., dipolar) for matrix elements to 
multiple quasiparticle-hole excitations, as we have seen 
it is for single quasiparticle-quasihole excitations, then 
these other contributions can be neglected. This seems 
likely to be correct, but as we do not have a proof, we 
will leave it as a conjecture that Eq. (4-6i) is an exact 
relation, which we call the "i<2 sum rule", and that it 
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holds for both the irreducible and reducible responses, as 
in generalized HF. If correct, we also obtain a relation of 
(1 + F^jm* to the LLL structure factor s(q) and V, as 



noted already in Sec. IV E 1 



When impurities are included, an improved approxi- 
mation is obtained by treating them diagrammatically 
similarly to the interaction lines as described at the be- 
ginning of this Section. In this Drude- (or SCBA-)-like 
approximation, the conductivity takes the same form as 
in Sec. [VE3. Based on the existing results [E3-p6|, 
we also expect that similar results hold for <fi > 1, with 
P=(2tt4,)-K 

We expect that the direct interaction of the particle 
with its correlation hole (or attached vortices) , described 
in Refs. [jl0| , p~2 24-2(| is contained in this description, but 
may not be easily obtainable diagrammatically. If it is 
obtained in some approximation, the effects stemming 
from F\ = — 1 will still be present when the approxima- 
tion is conserving. 

One other way that the FL picture could break down is 
by a pairing instability as in the theory of superconduc- 
tivity. The interaction in the quasiparticle-quasiparticle 
channel with quasiparticles of wavevectors k, — k can be 
considered using the ladder approximation. The dipolar 
nature of the coupling gives rise to an attractive interac- 
tion, as noticed by the authors of Refs. ]26| , |27|] . Since the 
system is compressible, this interaction is screened. In 
addition, the ladder series T, representing transverse and 
longitudinal gauge field fluctuations, can be exchanged 
between the fermions, and the transverse part can be 
combined with the interaction V . The transverse gauge 
field is believed to be pair-breaking when included in 
an Eliashberg equation treatment j fT6[ . The longitudi- 
nal part gives an extra repulsive short-range interaction, 
which also suppresses pairing, especially in the s-wave 
channel. Therefore the question of whether pairing is ac- 
tually expected to occur requires careful consideration. 
There is unpublished evidence that it does occur for 
bosons at v = 1 for some interactions p6| , ^7| . If pair- 
ing does occur, the system will become incompressible at 
low energies and long wavelengths, essentially because of 
the Meissner effect in the superfluid Fermi system: the 
diamagnetic susceptibility now behaves as Xd ~ — 1/9 2 > 
which inserted in our result for dn/ dp, shows the system 
is incompressible. This shows that it is not just the sym- 
metries of the Hamiltonian that make the ground state 
compressible in the FL-like state, but it is the fact that 
the state is assumed to be a normal (non-superfluid) liq- 
uid. 

Assuming the system is a FL, the scenario we have 
described here and in Sec. II D is essentially a FL cou- 
pled to an infinitely-strongly coupled gauge field (that 
represent F± = —1), with no CS term. The central point 
was the Ward identity that gave F\ = — 1. We con- 
nected this with the gauge invariance under XJ(N)r, or 
equivalently with conservation of G(q). Other authors 



have very recently commented on "translational invari- 
ance in momentum space" |2jj2^,^l|,^5[ , and its relation 
to some sort of gauge symmetry. We will try to make 
this more precise. The Hamiltonian Eq. (4.1) is invari- 
ant under shifts of the wavevectors of all the fermions by 
Q: k — > k + Q. The generator of a translation of the 
wavevectors of all fermions in such a system (or in ours) 
is 



d 2 k 



-kV k c k - (V k c k )c k 



(5.2) 



In first quantization and in position space, it is simply 
Y]. iv This is related to Galilean invariance in ordinary 
systems with finite bare mass tuq. If we rescale the gen- 
erator of Galilean transformation j77| to obtain shifts in 
ki instead of in v, = kj/mo, we obtain 



^(r 4 - tpi/mo), 



(5.3) 



and the second term can be dropped when mo — * oo. 
However, in this limit we obtain the ZFIM model, and the 
Galilean symmetry is enlarged to the local gauge symme- 
try generated by p(q), already discussed. In our system, 
by contrast, the gauge symmetry is generated by p fl (q), 



= N + 



-ijkAq t r 
2 e C k -iq C k4 



d 2 k 1 



(2?r)2 -q • [< A kc k c k 



4 v kCk - (V k c k )ck], 



(5.4) 



keeping all terms to linear order in q. Using the similar 
expansion of p L (q), the generator of shifts in k can be 
written as the first-order term in p R + p L . The other, un- 
used, pieces are the particle number N and the momen- 
tum Jkc^Ck, which are also conserved quantities (note 
that the terms in p L , p R linear in q are generators of mag- 
netic translations in the left, right coordinates, respec- 
tively, written in momentum space). Thus the "shifting" 
symmetry is part of the gauge symmetry, in combina- 
tion with other global symmetries, and not just part of 
the gauge symmetry as stated by SM. Even so, for some 
purposes, viewing it just as part of th e gauge symmetry 
can be useful, as we saw in Sec. [ID , and will again in 
the next paragraph. 

In unpublished work [ p7| , Haldane proposed to write 
the effective Hamiltonian of the quasiparticles, for the 
case of a finite system on a torus (say a square torus of 
side L), as 



1 



4m*N ^ 



(kj kj) 2 , 



(5.5) 



which possesses the shifting symmetry. In this system, 
shifting all the momenta by the smallest possible amount 
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2n/L changes the total momentum by 2irN/L, and gives 
a state equivalent to the original one 78 1 . The latter fact 
is assumed in numerical calculations, and such calcula- 
tions seem to confirm this form of the Hamiltonian. We 
may identify this Hamiltonian as similar to our 



A) 2 /(2m* 



(5.6) 



in the case of a spatially constant a + A, since (by an 
equation of motion) a + A = g/p — J^i kj/iV. The shift 
transformation is a gauge transformation (up to caveats 
just discussed) that does not change the physical states; 
this fact goes beyond the simple symmetry property pos- 
sessed by (5J5). Our Hamiltonian is preferable because, 
when a is allowed to vary spatially, it represents a local 
interaction, unlike Eq. fl5.5p . Integrating out a, and using 
the constraint on the density p R , we obtain a Hamilto- 
nian like that in SM, except that we have the effective 
mass m*, whereas in their work it appears at a stage 
where they instead have the bare mass m. This Hamilto- 
nian is also the starting point for the arguments of Ref. 



VI. CONCLUSION 

In this paper we have developed a truly lowest-Landau- 
level theory for the Fermi-liquid-like state of charged 
bosons at v = 1. We used a formalism of Pasquier and 



Haldane |52|, in which the composite fermion fields de- 
pend on two complex coordinates, one of which is the 
coordinate of the boson, and the other is in effect the 
coordinate of a vortex in the wavefunction of the other 
bosons, attached to the boson. The wavefunctions in 
both these coordinates are restricted to the lowest Lan- 
dau level, and there are operator constraints which fix the 
density in the vortex coordinates. The constraints imply 
that the system is a gauge theory. The effective theory 
for low-energy, long-wavelength phenomena is a Fermi 
liquid in which the fcrmions couple to a gauge field, for 
which there are no bare terms in the action. The ladder 



series treatment in Sec. [V, with the approximate form 
Eq. ( 4.44 ) , is equivalent to the RPA applied to this gauge 
field. Since there is no Chern-Simons term in the gauge 
field action, the longitudinal and transverse modes decou- 
ple. The longitudinal part, within RPA, gives rise to an 
effective scalar interaction at small momentum exchange 
of order the inverse density of states. This enforces the 
fixed-density constraint. The transverse part couples to 
the physical density, the first nontrivial term in which is 
dipolar in form and parity-violating. Each fermion car- 
ries a dipole moment equal to its wavevector. The result 
is a finite compressibility, and a low-frequency longitudi- 
nal conductivity that agrees with that in HLR. The gauge 
field obeys the same Chern-Simons equations relating it 



to the physical density and current as in the U(l) Chern- 
Simons fermion approach of HLR. Because there is no CS 
term in the action, the results nonetheless differ in form 
from those in the scenario for the fully-renormalized the- 
ory based on HLR. Although the gauge theory reduces to 
an ordinary U(l) theory at long wavelengths, this has to 
be supplemented by the expression for the density, which 
is a non-minimal coupling from the U(l) point of view. 
The form of the expression for the physical current in- 
timates that this is not the whole story, and we expect 
that the full gauge group will be involved in gen- 
eral. In view of existing results of other authors p4| , f25"| |, 
the results obtained here for <fr = 1 (bosons at v = 1) 
are expected to apply also for other cases of the FL-like 
state, when written in terms of p — (27T0) -1 and other 
parameters. There are many possible extensions and ap- 
plications of the present methods, to which we hope to 
return elsewhere. 
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APPENDIX A: NONCOMMUTATIVE 
GEOMETRY FOR PEDESTRIANS 

In this appendix we explain the formalism we use for 
states in and operators acting in the Hilbert space of a 
single particle in the lowest Landau level, in the simplest 
case of the infinite plane with uniform magnetic field, and 
magnetic length equal to 1 (see also Ref. @). This is 
equivalent to the "noncommutative plane" in noncommu- 
tative geometry. In particular we explain the "noncom- 
mutative Fourier transform" which we use extensively. 

The normalized basis states in coordinate representa- 
tion in the symmetric gauge are 

z m e -\\z\ 2 

u m (z) = , (Al) 
V ' V2^2™m! V ; 

A general state in the Hilbert space thus has wavefunc- 
tion ip(z) — /(z)e _ 4l z l , where / is a complex analytic 
function that does not grow too fast at infinity, so that 
/ \ip\ 2 is finite. All operators can be written as integral 
kernels, so that an operator a is represented by the kernel 
a(z, z'), which acts on states ip(z) as 



aip(z) = / d 2 z'a(z, z')tp(z'), 



(A2) 
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and matrix products become the "star product" a * b, the 
integral kernel of which is 



a * b(z, z') = / d 2 zi a(z, zi)b(zi, z r ). 



(A3) 

The operators themselves can, of course, be expanded as 

oo 

i,z')= ^2 a mn u m (z)u n (z'), (A4) 



alz, 



m.n—O 



so that a mn are elements of infinite matrices. 

Arbitrary operators in the larger Hilbert space of states 
in all Landau levels, that is all square-integrable complex 
functions in the plane (really, sections of the appropriate 
bundle), can be projected to the LLL. In particular, the 
identity 8(r— r') has matrix elements S mn in the orthonor- 
mal basis, and the corresponding operator as an integral 
kernel is 



S(z, z) = u m (z)u m (z') 



i-ex P (-i|z| 2 -i|zf + \zz'). (A5) 



As befits the identity, this obeys Sip = ip, 5*a = a*S = a. 
This operator also implements projection to the LLL. 

Another operator is defined by multiplication by the 
plane wave e lk r . Its projection to the LLL is 



d 2 z 1 <S(z,z 1 )e jk - ri <S(zi,z') 

S{z,z')e^ kz+k ^-^\ 



(A6) 



where, in this appendix, k = k x + ik y (elsewhere in the 
paper k = |k| for all vectors k). It is convenient to define 



r k (z,z') = 5{z,z')e^ {hz+ki,) -^ 2 . 



(A7) 



Thus fk = e , the adjoint of which is r_k, so 
Tk(z',z) = r_k(z,z'). The operator Tk has the effect of 
magnetic translation (i.e., translation which commutes 
with the Landau level index) by — ik or Ak in the plane 
J79[ . It obeys the well-known magnetic-translation rela- 
tion: 



Tk * T k ' = Tk+k'e 



(A8) 



Here \{kk' — kk') — ^ilmkk' = ijk A k', which is i 
times the (signed) area of the triangle formed by k, k', 
-(k + k'). 

The Tk are the natural functions for use in defining 
a "noncommutative Fourier transform" . The motivation 
is that functions (like the operator kernels) of z and z' 
are like wavefunctions for a single particle in zero mag- 
netic field, for which plane waves make sense. For such 
a function a(z, z'), we write 



a(z,z') = I ^a k T k {z,z'), 

Z7T 



and for the inverse transformation 



a k = / a * T_ k , 



(A9) 



(A10) 



where the integral is defined by J b — Tr b = J d 2 z b(z, z). 
The inversion theorem for this transform is easily proved 
by Gaussian integration. We note the orthonormality 
and completeness relations, 

J f k *f k , =27n5(k + k'), (All) 



cPk 



T k (z,z')T_ k (u;,w') = 5(z,w')5(w, z'). (A12) 



The "noncommutativity" of the transform shows up 
when one has convolutions, where the relation ( A8) must 
be used. 

In the main text the above formalism is applied to sec- 
ond quantized operators c, c', p L , p R , where it concerns 
their dependence on the z, w variables, and has nothing 
to do with the Fock space in which they act as operators. 
In the case studied in this paper, the Fourier transform 
can be applied to c and because the net magnetic field 
strength vanishes for v = 1/0 = 1. (For v ^ 1, one would 
require the full set of Landau- level states in the net, effec- 
tive magnetic field 0, projected to the z, w variables, in 
place of the plane waves which project to r k . The Fourier 
transform would still apply to p L and p R , of course.) For 
v = 1 wc define 



c(z, w) 



d 2 k 



-C k T k (z, to), 



(2tt) 3 / 2 
c k = (27T) 1 / 2 / c*f_ k ; 



(A13) 
(A14) 



the normalization has been chosen so as to obtain the 
conventional anticommutators in Eq. (3.14). For p L and 



p R we use the normalization given above for an arbi- 
trary a, and the properties of the r k 's lead to Eqs. ( p. 16 ), 
( |3.15 ). We also note that for the diagonal values z = z', 



p L (z,z) 



/(q) 



(2tt) 2 9 W 
< ><■" I d 2 rp L {z,z)e~ icir ', 



and similarly for p R . This exhibits the connection with 
GMP. 

Finally we note that other formulas of noncommutative 
geometry can be obtained in the integral kernel formal- 
ism. For example, the commutator in the star product, 



d*b~b*d= a* b , 



(A15) 
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defines the "Weyl-Moyal bracket" that generalizes the 
Poisson bracket of functions on the classical phase space 
to the quantum case. It is usually written as an infi- 
nite series of derivatives. Our integral kernel formula- 
tion avoids such series and allows generalization to other 
(e.g. compact) Riemann surfaces, or to nonuniform field 
strengths. In all cases, one can begin with an orthonor- 
mal set of LLL states, i.e. holomorphic sections of the 
appropriate bundle. A crucial operator is the "repro- 
ducing kernel" analogous to S(z,z'). This can be easily 
obtained for the sphere and the torus, for uniform field 
strength. 



APPENDIX B: HUBBARD-STRATONOVICH 
TRANSFORMATION AND THE l/M 
EXPANSION 

Here we show how to reproduce the results of the HF 
and ladder approximations as the saddle-point and Gaus- 
sian fluctuations in a Hubbard-Stratonovich field. First, 
one may replace the interaction term in the imaginary- 
time action by 



but not the Hartree, part of Hartee-Fock. Gaussian fluc- 
tuations in a around the saddle point reproduce the lad- 
der series. Thus the ladder series becomes the RPA in 
the a field. It should be possible to identify part of the 
a fluctuations as the gauge field, in a manner similar to 
that in some lattice models Q. 

In other problems, such a saddle point and Gaussian 
fluctuations are the leading terms in a l/M expansion, 
where M is the number of components of a field corre- 
sponding to our c, c 1 . We may introduce such compo- 
nents here, and then set M = 1 at the end, by replacing 
Cmn by c mna , where a — 1, . . . , M. The interaction is 
taken independent of a, so the system has SU(M) sym- 



f 4 

/ JJcPzi [c t (z 1 ,z 2 )c(z 3 ,z 4: )V(r2 - r 3 )a(z 4: , z 3 , z 2 , 



+ -\a(z 4 ,z 3 ,z 2 ,z 1 )\ 2 V(r 2 -r 3 )}, (Bl) 

(the r-dependence and r-integration is implicit) where 
a is a fourth-rank tensor field, written in the coordinate 
notation using LLL orthonormal functions as for c, , 
and is hermitian: 



cr(z 4 , z 3l z 2 , zi) = a(z 1 ,z 2l z 3 , z 4 ), 



(B2) 



and integrate functionally over a. Performing the latter 
functional integral reproduces the interaction term. The 
field a decouples the interaction in the exchange channel. 
The saddle point approximation for the a integral (along 
with the Lagrange multipliers) reproduces the exchange, 



metry. Then Eq. (Bl) now has the form 



M 

E 

a=l 



Y V(7- 



-M 



(B3) 



schematically. This appears suitable for l/M expansion, 
but there is a problem with the constraints. The latter 
must still be taken to be 



M 



EE 



(B4) 



in order to reproduce an M-component system of bosons, 
whatever the filling factor. To obtain zero net field for 
the fermions, we must be at total filling factor v = 1, 
so we must have p = l/2ir, that is of order M°, not M. 
Therefore not all the terms in the action are of order M, 
and we can expect problems with the l/M expansion. 
These are not necessarily completely fatal, however; an 
expansion can sometimes be obtained even in such cases 
(see Ref. |34j). It is not possible to rescale or redefine 
the model to avoid this problem. It could be avoided if 
we could attach l/M of a vortex to each particle (which 
would now be anyons, so that still creates fermions), 
as in the U(l) CS approach However, this is not 

possible in the present PH formalism. 
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